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Digital computer circuitry 


The possibility of making digital computers depends on 
the existence of non-linear devices capable in general of 
taking up two or more discrete states. Babbage’s pioneer 
attempts to make calculating machines in the nineteenth 
century,! for example, and indeed the design of all our 
modern desk calculators, depend essentially on the use of 
the discrete positions of toothed wheels to represent the 
digits 0 to 9 of the decimal system. Although a few early 
electronic computers also used the decimal system, and 
incorporated devices or groups of devices with the neces- 
sary ten stable states, it soon became apparent that greater 
speed and economy of operation could be achieved using 
the binary code of numbers, so that only two stable states 
needed to be used. The logical operations entailed in 
computation can also be conveniently specified by means 
of a binary switching algebra, so that exactly the same code 
can be used for instructions and for the numbers to which 
these instructions relate. 

All modern electronic digital computers, therefore, 
make use of binary physical devices. Their circuitry may 
be divided conveniently into two categories—a section 
which accomplishes the actual logical operations on 
numbers required by the calculation, and a section 
devoted to the storage of numbers and instructions at all 
the intermediate stages involved. Because the number of 
different operations is comparatively small, whereas the 
possible quantity of numerical data is almost unlimited, 
it has turned out that a fast and flexible computer needs a 
far larger number of binary devices in its store than it does 
in its logical circuits. To give an idea of the order of magni- 
tude involved, a large modern machine may store a 
million or more binary digits (bits) of information, while 
using only ten thousand or so binary elements in its logical 
circuitry. (Even this refers only to its rapid access store; 
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additional slow magnetic tape stores may be used to 
accommodate a hundred million or more bits.) Because of 
this, computers have in the past used quite different 
devices for the store and for the logical circuitry, a major 
consideration in the former being that of cost. In particular, 
whereas valve or transistor ‘flip-flops’ are widely used as 
high speed switches with inherent power gain in the logical 
circuits, their expense prohibits their use for the storage of 
information. This is now usually accomplished by means 
of magnetic ferrite cores, which although they are much 
slower in switching from one state to another, and do not 
possess power gain, are nevertheless at least an order of 
magnitude cheaper per bit. The presence or absence of 
available power gain in a device is often indicated by 
describing it as ‘active’ or ‘passive’, respectively. Passive 
devices, such as the ferrite core, must be used in conjunc- 
tion with separate amplifiers to boost up their output, 
whereas active binary devices are capable of switching one 
another independently. At the present time large ferrite 
core stores have been constructed, into and out of which 
information can be passed in roughly 1 microsecond, and 
it seems that this sort of speed is getting near the limit set 
by the physical properties of this device. Transistor flip- 
flop circuits for handling this information can be made 
which will switch at rather more than ten times this rate. 
Most computations, however, require a number of logical 
operations in between each reference to store; so that the 
relatively low speed of the latter is not such a handicap as 
might appear at first sight. 

Considerable efforts are being made in a number of 
laboratories to develop devices which will overcome the 
limitations of the present generation of computers, and it 
is the purpose of this review article to describe one such 
field of development, namely the utilization of the non- 
linear properties of superconductors in computer circuits. 
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Superconductivity 


Shortly after his successful liquefaction of helium at 
Leiden in 1908, Kamerlingh Onnes began a series of 
investigations? on the electrical resistance of metals in the 
new range of low temperatures which had thus been 
opened up. It was in the course of these experiments that 
in May, 1911 he first observed the abrupt drop in the resis- 
tivity of mercury at 4:15° K, and thus discovered the 
phenomenon of superconductivity. At the present day 
some twenty metallic elements and a large number of 
alloys are known to become superconductive, with tran- 
sition temperatures (T,) ranging from about 0-5° K (cad- 
mium, ruthenium, titanium) up to about 18° K for a 
compound of tin and niodium (Nb;Sn). 

The question as to whether the resistance of a metal in 
the superconducting state is truly zero or merely exceed- 
ingly small cannot yet be said to be finally settled. From 
the experimental point of view no trace of resistance has 
ever been measured, but of course this may mean only that 
the sensitivity of the experiments is not sufficiently refined ; 
the persistence of a current in the superconducting circuit 
known as the Crowe cell (see below) has proved that the 
maximum possible resistance is less than 4 x 107" of its 
normal value. The recent theories of superconductivity 
due to Bardeen, Cooper, and Schrieffer,> and Bogolyu- 
bov,* while remarkably successful in explaining many of 
the thermal and magnetic properties of superconductors, 
have so far failed to account satisfactorily for the dis- 
appearance of resistance (the phenomenon which gave 
superconductivity its name!), although they suggest the 
lines along which the search for the true explanation of 
this phenomenon should proceed. In the absence of such 
an explanation it is clearly impossible to say whether the 
theory of superconductivity will predict a genuinely zero 
resistance or not. 


Figure 1. H.-T curves of some superconductors 
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From the point of view of non-linear circuit elements 
the most important property of superconductors is the 
fact that their normal resistance can always be restored by 
the application of a large enough external magnetic field, 
which for a given pure superconductor is a well-defined 
function of temperature and is known as the critical field 
(H.). The critical field curves of a number of typical 
superconductors are shown in Figure 1, and the restora- 
tion of resistance in a thin mercury film as the applied 
magnetic field is increased in Figure 2. It is seen that, just 
as in the case of the change of resistance with temperature, 
this transition is an abrupt one. 
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Figure 2. Transition curves of mercury film 2,410 A thick 


For about twenty years after the discovery of super- 
conductivity it was assumed that the property of zero 
resistance, together with its restoration above the critical 
field, were sufficient to explain all the facts about super- 
conductors. Below H, any change in the applied field 
would induce currents in the surface of the metal just big 
enough to prevent any change of flux inside it, and because 
its resistivity was zero these currents would persist without 
any observable decay. The flux inside a superconductor 
would therefore remain constant at the value existing 
when the metal first went into the superconducting state. 
However, in 1933, Meissner and Ochsenfeld> showed con- 
clusively that if a superconductor was cooled in a magnetic 
field from above its transition temperature (so that the 
initial value of the flux inside it was not zero), then the onset 
of superconductivity was accompanied by a sudden expul- 
sion of flux from the interior of the specimen. This was 
shown by the ballistic galvanometer deflection observed 
as the specimen was cooled down, with a search coil 
surrounding it connected in series with the galvanometer. 
This classic experiment proved for the first time that it is 
not just the change in magnetic induction B which is zero 
in the interior of a superconductor, but the induction B 
itself. Instead of there being an infinite number of possible 
equilibrium magnetic states, only one, that with B=0, 
represents the true stable state of the material. A supercon- 
ductor thus behaves like a perfectly diamagnetic material 
of susceptibility « = —1/(4m), (since H+4nJ = B=0). 
This discovery of the Meissner effect, as it is called, 
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was of the greatest importance because it provided the jus- 
tification for the application of the theory of reversible 
thermodynamics to superconductors in a magnetic field, 
and enabled the correlation between their magnetic 
properties, such as critical field, and thermal properties, 
such as specific heat, to be worked out. 

An important corollary to the destruction of super- 
conductivity by an external magnetic field is its destruction 
by a current flowing in the superconductor itself. This was 
actually discovered by Onnes® before the discovery of the 
magnetic field effect, but is closely connected with it, the 
critical current in a superconducting wire of macroscopic 
dimensions being simply that current which will produce 
the critical field at the surface of the wire. The suggestion 
that this was the true explanation of the effect was first 
made by Silsbee’ and is known as the Silsbee hypothesis. 

One final property of superconductors needs to be 
mentioned before we pass on to a discussion of actual 
switching devices. The change in magnetic induction from 
zero inside a superconductor to a finite value outside, 
when the specimen is in an applied field, is not discontin- 
uous but takes place over a finite distance known as the 
penetration depth A, typically ~ 10~> cm. In other words, 
the surface currents, which keep the induction zero inside 
the metal, flow in a skin of finite thickness, and this is the 
penetration depth layer. Penetration effects become of 
particular importance in specimens in the form of thin 
metallic films whose thickness is of the same order of 
magnitude as the penetration depth. Two important 
changes then take place, an increase in the critical field 
required to destroy superconductivity in the film, and a 
decrease in its critical current. These effects will be des- 
cribed in greater detail when we come to consider devices 
based on superconducting thin films. 


Early superconducting switches 

The idea of using the superconducting transition as a 
switching mechanism is not a new one, and superconduct- 
ing switches have been used in laboratory apparatus for 
more than twenty years. One of the earliest of these was 
incorporated in the circuit of a superconducting galvano- 
meter designed by Grayson Smith and his coworkers at 
Toronto in 1935.8 This galvanometer consisted of a 
moving copper coil suspended in superconducting lead 
field coils, and was designed for the measurement of per- 
sistent currents induced in these field coils. The purpose 
of the superconducting switch was to introduce a small 
resistance in series with the field coils so that the persistent 
current dropped rapidly to zero, enabling the zero reading 
of the galvanometer to be obtained. The switch consisted 
simply of a short section of the superconducting circuit 
around which was wound a small copper coil, so that when 
a large enough current was passed through the latter the 
critical field of the superconductor was exceeded. A similar 
arrangement was used in conjunction with a much more 
sensitive superconducting galvanometer constructed by 
Pippard and Pullan,® except that their switch was an 
evaporated film of lead which could be heated above its 
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transition temperature by means of a heating strip of 
german silver. In this way a resistance of about 5 x 10-3 
could be introduced into the circuit, effectively ‘open 
circuiting’ a very small shunt across the superconducting 
galvanometer itself, and thereby enabling its sensitivity 
to be varied at will (a second resistance of about 2 x 10-°Q 
was permanently in series with the galvanometer). 

A further development of the superconducting switch 
was its use in superconducting ‘chopper’ amplifiers by 
Templeton!® and de Vroomen and van Baarle."! In these 
instruments the switch is driven alternately normal and 
superconducting by an alternating magnetic field so that 
a d.c. current in a circuit of very low resistance in series 
with the switch is modulated by it, producing an a.c. 
component which can be transformed up and amplified 
in a conventional way. The frequencies used by these 
authors were not high, 800 c/s and 10 c/s, respectively, but 
even so a serious difficulty was caused by pick-up due to 
stray fields from the coil. In spite of this both these ampli- 
fiers were capable of detecting signals as low as a few 
times.10-1)-V., 

The first deliberate attempt to apply the switching 
capabilities of superconductors to the development of 
computer components was made by Buck in 1956,!? who 
called the new device the ‘cryotron’. Like the switch used 
by Grayson Smith it depends on the destruction of super- 
conductivity by a magnetic field, and consists of a straight 
piece of superconducting wire, or ‘gate’, on which is 
wound a single layer coil through which a control current 
can be passed in order to drive the gate normal and thus 
restore its resistance. This control winding is itself made 
of a superconductor with a much higher transition tem- 
perature and critical field than the gate wire, so that even 
when it carries a current large enough to generate a field 
greater than the critical field of the gate it remains super- 
conducting itself. Suitable materials used by Buck for his 
cryotrons were tantalum for the gate (T, = 4-3° K) and 
niobium for the control winding (7, = 8° K); at the normal 
boiling point of liquid helium (4-215° K) the critical field 
of the gate was then about 40 oersteds. 

An important feature of the cryotron is that, by a 
suitable choice of the dimensions of the gate and control 
circuits, it is easy to arrange for the critical current i, of 
the gate wire to be larger than the current i, in the control 
winding needed to restore the gate resistance, or ‘critical 
control current’. Then if H, is the critical field of the gate 
wire we have (by Silsbee’s hypothesis) 


H, = 4i,/d be a) 
and H, = 4nni, dps 4) 


where dis the diameter of the gate wire and n is the number 
of turns per unit length of control winding. Thus 


ili, = md = K ve) 


The number K may be called the ‘current gain’ of the 
cryotron. For the cryotron dimensions used by Buck (gate 
diameter 0-009 in., control winding 250 turns per inch of 
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0-003 in. insulated niobium wire), K =7. In practice the 
observed value of current gain was usually somewhat 
smaller than that calculated, presumably due to small flaws 
in the tantalum gate wire. Provided that K is greater than 
unity it is possible to construct switching circuits with 
identical cryotrons (using superconducting interconnec- 
tions) without any additional power amplification, in 


Schematic 


Figure 3. Cryotron flip-flop 


just the same way as is done with relays. The value of K also 
fixes the allowable tolerance on the cryotron parameters, 
for a circuit in which there is a statistical spread of these 
parameters will always work satisfactorily provided that 
the ratio of the highest to the lowest critical current is less 
than the smallest current gain. Since in the steady state 
the current in such a circuit will always flow along a 
completely superconducting path there is no power dissi- 
pation and heat is generated only during the actual switch- 
ing process. The basic simplicity and universality of the 
cryotron for the purpose of switching circuitry are almost 
unique among digital devices, and it is therefore not 
surprising that they have, as we shall see, attracted so much 
attention in recent years. 

The simplest example of a cryotron switching circuit is 
the bistable element or flip-flop, shown in Figure 3, con- 
sisting of two cryotrons A and B connected in parallel with 
the gate of one in series with the control circuit of the other. 
The current 7 supplied externally can flow either through 
gate A and control winding B, holding gate B resistive, or 
through gate B and control winding A, holding gate A 
resistive. It is, of course, necessary that i, > 7 > i,. Further- 
more this circuit as it stands provides no means of con- 
trolling the current path: two further cryotrons C and D 
(Figure 4) are needed to act as ‘write-in’ triggers, and one 
more E as a ‘read-out’ circuit. The complete flip-flop 
therefore consists of five cryotrons. 


Figure 4. Complete flip-flop 
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The rate at which the current i can be switched from one 
path to the other depends on the self-inductance of the 
whole loop and on the normal gate resistance R associated 
with each cryotron. If most of the self-inductance is due 
to the control winding, as is usually the case, it turns out 
that the overall switching time is about 7L/R, where L is 
the self-inductance of one control winding, and in the case 
of the cryotrons described above this time is of the order 
of 100 psec. A further limitation to the speed of the device 
is set by the heat generated during switching while current 
is still flowing momentarily through a resistive arm of the 
flip-flop. Since the rate of heating is Ri? and since the time 
for which it continues is of the order of L/R the total heat 
generated is roughly Li’. This heat warms up the gate wire, 
and it is found in practice that, because the maximum 
permissible temperature rise of the tantalum is less than 
0-1° K, this effect limits the frequency at which the flip-flop 
may be switched to about 1,000 c/s. 


Lead stri 
(control 


Silicon monoxide 
insulation 


Figure 5. Strip cryotron 


We see, therefore, that the Buck wire-wound cryotron is 
severely limited in its speed of operation, mainly by its 
comparatively large self-inductance. The efforts made in 
the last few years to overcome this limitation, and the 
extent to which they have been successful will be described 
below. 


Fast superconductive devices 


A simple modification of the basic design of the cryo- 
tron in order to reduce its self-inductance, and at the same 
time to increase its normal resistance, immediately suggests 
itself. This is to make both gate and control circuits of thin 
evaporated superconducting films, instead of macroscopic 
wires, in the way shown schematically in Figure 5. In this 
construction a current through the lead control strip gives 
rise to a magnetic field between it and the ground plane, 
so that the narrow region of the tin gate lying immediately 
beneath the control strip can thereby be driven normal. 
The geometrical design of this device has given rise to the 
name ‘planar cryotron’ to distinguish it from Buck’s 
original invention. It was first tested at the laboratories of 
1.B.M. in America, but has shown so many promising 
features that many other industrial laboratories in the 
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U.S.A. have started working vigorously on it, though it has 
so far been almost entirely neglected by industry in the 
U.K. 

The chief problem encountered in fabricating the 
planar cryotron has proved to be that of achieving a 
useful current gain. The first tests were made on twocrossed 
strips with no ground plane, and showed not only that the 
gate in this arrangement had a disappointingly low critical 
current but also, from the width of its transition as the 
control current was-aised, that the magnetic field due to 
the latter was too spread out. The presence of the ground 
plane not only concentrates this field (and incidentally 
thereby considerably reduces the self-inductance of the 
circuit), but probably also makes the current distribution 
in the gate much more uniform across its width, thereby 
raising its critical current. If the widths of the gate and 
control strips are w, and w,, respectively, then the fields 
beneath them produced by unit current are 47/w, and 
47/w, and, assuming again that Silsbee’s hypothesis holds 
in this case, the current gain K of the cryotron should be 
W,/W,, or width ratio W. In practice it is found extremely 
difficult even using a ground plane to achieve this ideal 
current gain, so it is useful to introduce a cryotron 
efficiency e defined by e = K/W. 

The I.B.M. workers have so far been the most successful 
in attaining values of e close to unity, and have traced one 
particular cause of trouble which is of some interest. The 
gate strip is usually made by evaporating tin (or indium) 
on to the substrate (after the ground plane and silicon mon- 
oxide insulation), through a suitable mask, and indoing so 
the edges of the strip always show a ‘penumbral’ region in 
which the thickness of the strip gradually falls off to zero. 
Now, as was mentioned above, such a region of very thin 
film has a critical field greater than the bulk metal, so that 
it will remain superconducting, and thus short out the 
central region of film, after the latter has been driven 
normal by the control current. The presence of the 
penumbra thus increases the effective critical control 
current and thereby decreases the current gain, and it 
appears that its removal, either by physically scoring 
through it or by other means such as depositing on a 
pre-heated substrate (which has the effect of causing the 
penumbra to disintegrate into isolated ‘islands’ of film), 
can increase the current gain by a factor of 2 or 3. 

The increased critical field of a very thin film is, as we 
have seen, due to the penetration of the magnetic field 
into the surface layer of the superconductor, and it sets a 
lower limit to the useful thickness of the gate in the planar 
cryotron. If current gain were the only important factor 
there would be no point in making the films as thin as 
possible, but considerations of switching speed show that 
this is not so. Once again it is L/R which is the important 
factor. If the main contribution to self-inductance in a 
cryotron circuit is that of the control strip, then if the 
insulation thickness is d cm and the gate thickness and 
resistivity are t¢cm and p QQcm, respectively, we have 


L = 47x 10~° dw,/w, H - «+ (4) 
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and R = pw,/tw, Q asd5) 


so that L|R = 47 x 107? dtW?/p sec 


= 47 x 10~? dtK?/pe? sec ... (6) 


For a given current gain L/R has a minimum value for 
some thickness ¢ usually somewhat greater than the pene- 
tration depth of the superconductor, owing to the rapid 
decrease in e as t is reduced. The theory of the destruction 
of superconductivity in thin films both by external fields 
and by currents is not yet sufficiently well understood to 
say whether the best observed values of e are ‘ideal’ or 
still limited by imperfections in the film. However, for 
gate thicknesses of about 8,000 A, efficiencies of over 0-9 
have been reported by Rosenberger,'? so that no order of 
magnitude improvement can be expected. The value of 
L/R predicted by equation (6) is then about 2 x 10~® sec 
for K = 4, p= 10~7 Qcm, d= 1 micron, and agrees quite 
well with the observed switching time (7L/R) of an actual 
flip-flop. 

Further improvements in the speed of the planar cryo- 
tron may be expected from developments in two directions. 
First one may expect a significant improvement in the 
figures quoted above for the factors in equation (6)—the 
use of alloys may enable p to be increased, the insulation 
thickness may be further reduced, and it may prove pos- 
sible to decrease ¢ without a corresponding deterioration of 
e. Better control of evaporating conditions will probably 
result in cryotrons of greater reproducibility, thus enabling 
a smaller current gain to be used. Conventional improve- 
ments of this kind will probably resultina reduction of L/R 
by a further factor of 5 to 10. In the second place a recent 
development in the geometrical design of the planar cryo- 
tron at I.B.M. may lead to even higher speeds. This is the 
use of what I.B.M. call the ‘in-line’ cryotron, in which the 
control and gate strips, instead of crossing one another at 
right angles, are actually parallel, so that the entire length 
of gate lying below the control strip is driven normal, 
instead of just a narrow section. At first sight this geometry 
would seem inapplicable, since the current gain must 
inevitably be less than unity. However, by the use of a 
further biasing control strip placed above the ordinary 
control, advantage can be taken of the fact that with 
parallel gate and control their magnetic fields can be either 
parallel (currents flowing in the same direction) or anti- 
parallel (currents opposed). A steady biasing field can thus 
be made to increase the control field but decrease the self- 
field due to the gate, and a current gain greater than unity 
can thus be achieved. The in-line cryotron has a normal 
gate resistance very much larger than the conventional 
design, and it is hoped ultimately that it can actually be 
matched approximately to the characteristic impedance 
of the circuitry linking the cryotrons together. This cir- 
cuitry consists of evaporated superconducting lead strips 
just like the control strips of the cryotrons themselves, and 
forms, with the ground plane, a set of strip transmission 
lines of roughly | Q characteristic impedance. At present 
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it requires an in-line cryotron gate about 4 cm long to 
match this resistance, but this length can probably be 
considerably reduced. If such matching can be achieved 
the speed of switching of cryotron circuits will be increased 
to such an extent that it becomes limited mainly by the 
propagation time along the circuit at the velocity of 
electromagnetic waves in the dielectric insulation. 

The planar cryotron is not the only fast superconductive 
computer device to appear in the last few years. It was 
mentioned earlier that one of the limitations of present 
computers is the comparatively slow rate of access to the 
store set by the switching speed of the ferrite core. In the 
quest for higher speed, therefore, considerable attention 
has been devoted to the development of faster passive 
devices suitable for information storage. One such device 
which shows promise makes use of the possibility of setting 
up persistent currents in a superconducting circuit, the 
direction of these currents representing the two states 
required to store information. It was invented by Crowe'* 
at I.B.M. and has now been studied in some detail both in 
the U.S.A. and also in the U.K. at R.R.E. and at S.E.R.L. 

The Crowe cell consists of a thin film of tin or lead of the 
shape shown in Figure 6(a) on a glass or mica substrate. 
The two states of the device are represented by currents 
flowing in the directions shown in Figure 6(b) and 6(c), 
and are set up by current pulses of opposite polarity and 
appropriate amplitude in the drive wire, which is closely 
linked inductively with the superconducting circuit. 
Changes of state are detected by means of the pick-up wire 
beneath the ‘crossbar’. Because the superconducting 
circuit acts as a perfect screen for flux through the two 
D-shaped holes, a signal is only obtained when the cross- 
bar is driven normal, and it can be shown that this occurs 
only when the drive current pulse is of opposite polarity 
to the previous one. Thus the ‘information’ which is 
stored in the cell can always be read out destructively, in 


Drive wire 


x 
Evaporated 
metal film 


Figure 6. Diagram of Crowe cell: (a) General view; (b) ‘0’ state; 
(c) ‘1’ state 
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the same way as with a ferrite core, by ‘interrogating’ it 
with a current pulse of the appropriate polarity to drive it 
into the ‘0’ state. A previously stored ‘1’ will then give an 
output signal, but a previously stored ‘0’ none. 

The Crowe cell depends on the destruction of super- 
conductivity in the crossbar when the induced current in 
it exceeds the critical value. It differs in this respect from 
the cryotron whose operation depends upon the external 
application of the critical magnetic field of the gate. For 
this reason there is no inherent disadvantage in using 
much thinner films for the Crowe cell, and the switching 
speed can be thereby increased. As usual, this depends 
mainly on the L/R time constant of the circuit when the 
normal resistance of the crossbar is restored; and the self- 
inductance is in fact appreciably larger than in the planar 
cryotron, about 10~? H. However, the crossbar resistance 
can be quite large, a few tenths of an ohm, so that the 
switching speed of the Crowe cell is high. The best figure 
reported so far is about 3 x 10~? sec, or some 100 times 
faster than the highest speed attainable with the ferrite 
core. With this high speed, with its fairly small drive 
current (~ 100 mA), and with the possibility of evaporat- 
ing it in large arrays at low cost, the Crowe cell offers 
attractive possibilities for a large rapid access computer 
store. Whether these can be fulfilled or not depends largely 
on the question of reproducibility, for the tolerance allow- 
able on the critical current of the crossbar is much smaller 
than is the case for the gate of the planar cryotron. 


Cryotron circuitry 


Those who have worked on the planar cryotron are so 
confident of its ultimate success that a considerable amount 
of work has been devoted to the design of cryotron 
computer circuits, well before the cryotrons themselves 
have become available in appropriate numbers. It is 
important to stress here that the ultimate goal in this field 
is the construction of a computer whose circuits, apart 
from the input and output, are built exclusively with 
cryotrons; no other type of circuit element is either needed 
or desired. 

It has been pointed out by Haynes!® that the cryotron 
has the unique property among digital computer devices 
of providing a means of storing information without power 
dissipation, and at the same time of yielding a d.c. indi- 
cation of its state. This can be achieved, as we shall see, 
by using the cryotron in a persistent current loop, as in 
the Crowe cell, but at the same time causing the persistent 
current to act as the control for a second cryotron. All 
other types of fast storage device belong to one of two 
categories: either they give a d.c. indication of their state, 
but require power dissipation to maintain that state (for 
example the transistor flip-flop), or they require no 
continuous power dissipation but yield no d.c. indication 
of their state (for example, the ferrite core). 

A single cryotron performs exactly the same logical 
function in a circuit as a simple relay with a normally 
closed contact. Thus, when no current flows in the control 
strip, then current may flow in the gate and vice versa. 
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Just as in the case of relays, any logical function whatever 
may be built up by a suitable combination of cryotrons. 
This similarity between cryotrons and relays enables 
computer circuitry to be designed on an asynchronous 
principle, or in other words without the need for a syn- 
chronizing clock pulse generator. One can set up a network 


Figure 7. Cryotron persistent current store 


of cryotron gates by means of the control currents appro- 
priate to any chosen input function, and then apply a 
probing current to the entire network of cryotrons. This 
probing current will distribute itself through the network 
as rapidly as possible in a pattern of paths which is 
switched successively from one arrangement to another 
until it attains the stable state yielding the corresponding 
output function. 

It is clearly of primary importance to devise a 
cryotron circuit suitable for a large random access 
store, provided with a convenient method of ad- 
dress selection. A solution to this problem has 
been given by Haynes!® by means of the circuit 
of Figure 7 which requires three cryotrons per bit 
of storage, and uses the presence or absence of a 
persistent current in the loop abcd to represent a 
‘1’ ora ‘0’, respectively. Figure 7 shows a portion 
containing 4 bits, of a two-dimensional store 
arranged in rows of words so that a given column 
contains bits of the same degree of significance in 
all words. Two-dimensional selection for writing 
and completely separate non-destructive reading 
is provided by ‘input’ and ‘write select’, ‘sense’ 
and ‘read select’ lines, respectively. 

To storea ‘1’ inthe loop abcd, current is passed 
through the appropriate digit input line D and 
the write select line E for word 1 is pulsed. This 
drives the gate of cryotron B resistive so that the 
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input current must switch to the shunting loop abcd, of 
appreciably higher self-inductance than the direct path 
ad. Flux is now linked with the loop abcd, and on remov- 
ing the currents from E and then D a persistent current 
will remain in the loop, large enough to control the gate 
of cryotron C. To store a ‘0’ or to switch a ‘1’ to a ‘0’ 
the write select line is pulsed with no current flowing in 
the input line. 

Similarly, to read the contents of loop abcd current is 
applied to the digit sense line F while the read select line 
G is pulsed. The latter drives the gate of cryotron A 
resistive, and, if the loop abcd contains no persistent 
current so that gate C is superconducting, the sense current 
is switched to the shunting loop hgfe but continues to flow 
in line F. If, however, a ‘1’ is stored in loop abcd, gates 
A and C are both resistive and the current in F is switched 
to some other superconducting path in parallel with F (not 
shown in Figure 7), thereby controlling a digit register 
cryotron which will then sense the presence of the stored 
“1”. The separation of the read and write selection circuits 
in this way enables simultaneous reading of one word and 
writing of another to be performed if necessary. The same 
basic storage unit, in conjunction with slightly more 
complicated access circuits, makes it possible to perform 
more sophisticated logical operations such as transfer of 
data from one word or register to another, right shifting in 
storage, or interchange of data between two registers. 

As a final example of cryotron circuitry we give the 
binary full adder for a single bit.'® This circuit requires 
fourteen cryotrons, and is shown in Figure 8. A full adder 
has three inputs, namely the carry from the previous stage 
and the two digits to be added, and two outputs, the sum 
and the carry to the next stage. These are indicated in 
Figure 8 by the letters Co, A, B (the inputs), S;,and C, (the 
outputs), respectively. A bar over the letter indicates a ‘0’ 
and its absence a ‘1’. The first stage of the operation 


S 
'B Intermediate result from (a) 
fol 101 [10 
ofotr a Beata 
i] 1 [10] | Jor fio} | 


(a) (b) 
Figure 8. Binary full adder 
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involves acting on the input carry by the B input and must 
produce the intermediate sum shown in Figure 8(a) which 
has the three possible values, 0, 1, and 10. Thus the two 
lines corresponding to Cy. B and Cy. B (remembering that 
the presence of control current in a cryotron denies 
current to its gate line) are recombined, giving three inter- 
mediate result lines. These must now split into six so that 
either A or A can operate to produce the required full sum 
shown in Figure 8(b), in which the first digit is the C; carry 
and the second the S; sum. Again recombination takes 
place to give the four possible results 00, 01, 10, and 11 on 
four lines which generate the sum output via the gates in 
lines S; and S, and themselves recombine into the two 
carry lines C, and C;. 


The future of cryotrons 


Of all the factors which determine the success or failure 
of a new invention such as the cryotron the most important 
is probably also the most difficult to assess—namely the 
economic factor. As regards its technical capabilities the 
planar cryotron seems to be ideally suited for future 
mathematical computers; it combines simplicity, versa- 
tility, compactness, high speed, negligible power dissi- 
pation, and, above all, ease for translating the logical 
requirements of the computer designer into an actual 
circuit, since it is a digital device par excellence, needing 
no other ancillary circuit elements in combination with it. 
The need for a liquid helium refrigerator in a cryotron 
computer should be no barrier, since modern cryogenic 
techniques have shown that such machines can be made 
reliable over running periods of over 1,000 hr.!’ If it were 
economically feasible to spend as much money on the 
development of the cryotron as has been spent on, say, 
the transistor it might lead the field in computer design. 
But, of course, the economics of the transistor depend, not 
oncomputers, but onits tremendous market incommercial 
‘radio and television, among other things, and it is this 
which has enabled computer designers to cash in on the 
availability of comparatively cheap high speed transistors. 
The demand for large fast mathematical computers, 
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particularly in the U.K., is a very restricted one, and 
unfortunately some of the chief customers, the univer- 
sities, are not able to buy as many as they would like 
because of lack of funds. The electronics industry is, in 
fact, loth to sink capital on a completely new project 
such as the cryotron, restricted entirely to the computer 
field. 

Nevertheless, when one considers the enormous number 
of possible applications in an industrial society of com- 
puters, data processing machines, and automation in 
general, it would seem economically unrealistic to neglect 
a device which already, in its infancy, shows so much 
promise, and which could find almost universal appli- 
cation in machines of this kind. 
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Abstracts 


THE VARIATION OF THE DENSITY OF LIQUID 
NITROGEN AND LIQUID OXYGEN AS A 
FUNCTION OF PRESSURE 


A. VAN ITTERBEEK O. VERBEKE 


MEASUREMENTS have been carried out on the variation of the specific 
gravity of liquid nitrogen .and liquid oxygen up to 850 kg/cm?. The 
results can be represented analytically by means of the following 
equation (p in g/cm’, P in kg/cm?): 
N,:90:26°K; p = 0:7497+ 3-493 .10-4P— 3-431.10-7P2+ 
+0-01713.10-8P? 
77:31° K; p = 0:8079+ 2-521 .10-4P— 2-031. 10-7P?+ 
+0-007995 .10-8P® 
O,:90:07° K; p = 1:1419+ 1-969. 10-4P— 1-:233P?.10-7+ 
+ 0:4682.10-1°P 
77:35°K; p = 1:1997+ 1-608. 10-4P—0-8078P?. 10-7 + 
+ 0-298 . 10-19P3 


THE DENSITIES OF SATURATED LIQUID 
HYDROGEN 


D. E. DILLER 
L. A. WEBER 


R. D. GOODWIN 
H. M. RODER 


Densities of para-hydrogen from 17° to 32° K are obtained by 
extrapolating newly determined compressed liquid isotherms to 
corresponding vapour pressures. Results are combined with earlier 
data for the range 14° to 20°K and compared with normal hydrogen 
by means of density-temperature relationships. 


SMALL SCALE HYDROGEN LIQUEFACTION 


B. YATES F. E. HOARE 


A DESCRIPTION is given of the procedure adopted in designing a 
hydrogen liquefier to produce a required output. The procedure is 
applied to a case where 5 |. of liquid hydrogen were required to be 
produced without the use of a compressor. Constructional details of 
the purification and liquefaction systems of an accessible and easily 
constructed machine are included, together with performance 
figures. 


TOLERANCES ASSOCIATED WITH THE 
OPERATION OF THE CROWE STORAGE CELL IN 
THE COINCIDENT CURRENT MODE 


R. V. PEACOCK 


Tue coincident current mode of operation of the Crowe storage cell is 
considered. The relationship between the tolerances on the device 
parameters and the drive currents is calculated as a function of the 
shape of the adiabatic curve of resistance against current of the 
crossbar. 
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A METHOD FOR MEASURING THE HEAT 
CAPACITY OF CONDENSED GASES ABOVE 
THEIR BOILING POINT 


A. V. VORONEL P. G. STRELKOV 


A METHOD has been developed for measuring the heat capacity with 
stirring at pressures up to several tens of atmospheres. The heat 
capacity of the calorimeter, which is approximately 105 cm? in 
volume and can stand 70 atm, amounted to 13-20 per cent of the 
measured quantity. The magnetic stirrer working at low temperature 
and high pressure, does not introduce any noticeable heating and 
reduces the time for achieving equilibrium by 1-5 to 2 times. Trial 
measurements were made of the heat capacity of nitrogen in the 
critical region. 


THE MEASUREMENT OF THERMAL 
EXPANSION OF METALS AT LOW 
TEMPERATURES 


K. ANDRES 


AN apparatus for measuring thermal expansion at low temperatures 
is described. By means of an optical technique, it is possible to 
detect length changes of 1 A (10-8cm) in a specimen 7 cm long. 
Results are given for the lattice and electronic contributions to 
thermal expansion in the metals Al, Pb, Ta, Nb, Mg, and W. 


CARBON RESISTANCE THERMOMETERS 
FOR LOW TEMPERATURES 


N. N. MIKHAILOV A. YA. KAGANOVSKII 


THE technique of preparation from petroleum pitch and the pro- 
perties of carbon thermometers, for measuring low temperatures 
with high sensitivity and satisfactory reproducibility are described. 
The form of the temperature dependence of the resistance of these 
thermometers is determined by the firing temperature. 


A HYDROGEN LIQUEFIER WITH AN OUTPUT 
OF 50 1./hr OF LIQUID HYDROGEN 


A. G. ZEL’DOVICH YU. K. PILIPENKO 


THE principal of operation and construction of a liquefier producing 
50 1./hr of normal, or 30 1./hr para-hydrogen is described. The 
possibility is noted of the operation of the liquefier in a refrigerating 
cycle. A table showing the performance of the liquefier is given, and 
also the results of measuring the losses. 


THE EFFECT OF THE MARTENSITIC 
TRANSFORMATION ON THE ELECTRICAL 
RESISTANCE OF LITHIUM AND DILUTE 
LITHIUM-MAGNESIUM ALLOYS 


J. S. DUGDALE D. GUGAN 


THE effect of the martensitic transformation on the electrical re- 
sistivity of lithium and dilute lithium—magnesium alloys has been 
studied. The resistivity of pure lithium at 80° K decreases by approxi- 
mately 0:5 per cent on transforming, and that of the alloys by a 
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great deal more. By means of experiments using thermal cycling and 
mechanical deformation we find that the intrinsic residual resistivities 
of the h.c.p. and f.c.c. phases of lithium and the magnesium-lithium 
alloys are, respectively, about 10 and 15 per cent greater than that of 
the b.c.c. phase. The experiments show that there are large departures 
from Matthiesen’s rule in all three phases of the alloys and in at least 
the b.c.c. phase of quite pure lithium. These departures account for 
the much larger changes in resistance that are observed when the 
alloys transform than are observed when lithium transforms. 


Resumes 
LA VARIATION DE LA DENSITE DE 


L’AZOTE ET DE L’OXYGENE LIQUIDE EN 
FONCTION DE LA PRESSION 


A. VAN ITTERBEEK O. VERBEKE 


ON a fait des mesures de variation de masse spécifique d’azote liquide 
et d’oxygéne liquide jusqu’a 850 kg/cm?. Les résultats peuvent étre 
représentés a l’aide des équations suivantes (peng/cm’, P enkg/cm?): 


Nz: 90:26°K; p = 0:7497+3-493. 10-4 P—3-431.10-7 P?+ 
+0-01713.10-® P® 
77-31°K; p = 0°8079+2:521.10-* P—2-031. 10-7 P?+ 
+0-007995.10-§ P® 
O2: 90:07°K 3p = 1-1419+- 1-969. 10-4 P— 1-233, P?. 107+ 
+0-4682.10-1° PS 
77:35°K; p= 1-1997+1-608 . 10-* P—0:8078 P?. 1077+ 
0298.) 10mUP* 


LA DENSITE DE L’HYDROGENE LIQUIDE 
SOUS SA PRESSION DE VAPEUR SATURANTE 


R. D. GOODWIN D. E. DILLER 
H. M. RODER L. A. WEBER 


Les densités d’hydrogéne para entre 17° et 32°K ont été obtenues en 
extrapolant aux pressions de vapeur correspondantes les isothermes 
du liquide comprimé nouvellement déterminées. Les résultats sont 
combinés avec des valeurs antérieures dans l’intervalle de 14 4 20°K 
et comparés avec I’hydrogéne normal a l’aide des relations densité- 
température. 


LIQUEFACTION DE L’HYDROGENE A 
PETTIE ECHELCE 


B. YATES F. E. HOARE 


DESCRIPTION du procédé adopté pour réaliser un liquéfacteur de 
débit donné. La méthode est appliquée ici, pour obtenir 5 litres 
d’hydrogéne liquide sans l’aide de compresseur. On trouvera dans 
ce texte, les détails de construction des systémes de purification et 
de liquéfaction d’un appareil accessible et facile 4 construire, ainsi 
que ses performances. 


TOLERANCES ASSOCIES AVEC L’OPERATION 
DE L7-ELEMENT D’EMMAGASINAGE CROWE 
DANS LE MODE PAR COURANT 
COINCIDENT 


R. V. PEACOCK 


Le mode d’opération par courant coincident de |’élément d’em- 
magasinage Crowe est considéré. Le rapport entre les tolérances 
sur les paramétres de l’appareil et les courants de transmission 
est calculé comme fonction de la forme de la courbe de résistance 
adiabatique contre le courant de la traverse. 
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A LIQUID HELIUM LEVEL INDICATOR 
A. E. ROVINSKII 


Tue idea of using the phenomenon of the transition of pure tantalum 
into the superconducting state at 4-3°K for measuring the level of 
liquid helium was proposed by Feldmeier, and first applied by Rasor 
in an apparatus for continuous measurement of the quantity of 
liquid in a closed vessel. In a liquid level indicator with a feeler in the 
form of a probe it is desirable to reduce the dimensions and mass of 
the immersed tantalum element, so that it is advantageous to change 
the electrical circuit of the instrument somewhat in order to maintain 
a sufficiently high sensitivity. 


UNE METHODE DE MESURE DE LA 
CAPACITE THERMIQUE DE GAZ CONDENSES 
AU DESSUS DE LEUR POINT D’EBULLITION 


A. V. VORONEL P. G. STRELKOV 


ON a mis au point une méthode de mesure de la capacité thermique, 
avec agitation, a des pressions allant jusqu’a plusieurs dizaines 
d’atmosphéres. La capacité thermique du calorimétre, d’un volume 
d’environ 105 cm? et pouvant résister 4 70 atmosphéres, était de 
13-20 pour cent de la quantité mesurée. L’agitateur magnétique, 
fonctionnant a basse température et haute pression, n’a pas créé de 
chaleur appréciable et a réduit le temps nécessaire pour obtenir 
l’équilibre de 1,5 a 2 fois. Des mesures expérimentales de la capacité 
thermique de l’azote dans la région critique ont été réalisées. 


MESURE DE LA DILATATION THERMIQUE 
DES METAUX A BASSE TEMPERATURE 


K. ANDRES 


DESCRIPTION d’un appareil permettant de mesurer les dilatations 
thermiques a4 basse température. Une technique optique permet de 
détecter des variations de longueur de 1 A (10-* cm) sur un échan- 
tillon de 7 cm de long. Les contributions a la dilatation thermique 


des électrons et du réseau sont données pour les métaux Al, Pb, Ta, 
Nb, Mg, et W. 


THERMOMETRE A RESISTANCE DE 
CARBONE POUR BASSES TEMPERATURES 


N. N. MIKHAILOV A. YA. KAGANOVSKII 


DESCRIPTION d’une technique de préparation a partir des goudrons 
de pétrole et propriétés des thermométres a carbone pour la mesure 
avec grande sensibilité et reproductibilité satisfaisante des basses 
températures. La variation de la résistance de ces thermomeéetres en 
fonction de la température dépend de la température de combustion. 


LIQUEFACTEUR D’HYDROGENE DE DEBIT 
50 1./hr EN HYDROGENE LIQUIDE 


A. G. ZEL’DOVICH YU. K. PILIPENKO 


Les auteurs décrivent l’essentiel de l’utilisation et de la construction 
d’un liquéfacteur produisant 50 1./hr @hydrogéne normal et 30 
l./hr de parahydrogéne. Le liquéfacteur peut également travailler 
en cycle réfrigérant. Un tableau montre les possibilités de ce 
liquéfacteur et le résultat de la mesure des pertes. 
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EFFET DE LA TRANSFORMATION MARTEN- 
SITIQUE SUR LA RESISTANCE ELECTRIQUE DU 
LITHIUM ET DES ALLIAGES DILUES 
LITHIUM-MAGNESIUM 


J. S. DUGDALE D. GUGAN 


Nous avons étudié l’effet de la transformation martensitique sur la 
résistivité électrique du lithium et des alliages dilués lithium—mag- 
nésium. La résistivité du lithium pur, a 80° K, décroit d’environ 0,5 % 
lors de la transformation, celle des alliages dilués diminue encore bien 
plus. Des expériences de cyclage thermique et de déformation mécan- 
ique nous ont permis de trouver que les résistivités intrinséques 
résiduelles des phases hexagonale compacte et cubique faces centrées 
du lithium et des alliages lithium—magnésium sont, respectivement, 
10 et 15% plus élevées que celle de la phase cubique centrée. Les 
expériences montrent, pour les trois phases des alliages et au moins 
pour la phase cubique centrée du lithium trés pur, de larges écarts 
a la régle de Matthiesen. Ces écarts expliquent que les variations 


Auszuige 


SCHWANKUNGEN DER DICHTE VON 
FLUSSIGEM STICKSTOFF UND SAUERSTOFF 
ALS FUNKTION DES DRUCKS 


A. VAN ITTERBEEK O. VERBEKE 


Der Artikel bespricht Messungen der Dichte von fliissigem Stick- 
stoff und Sauerstoff bei Driicken bis zu 850 kg/cm?. Die Resultate 
k6nnen analytisch durch die folgenden Gleichungen wiedergegeben 
werden (p in g/cm, P in kg/cm?): 
N,2: 90:26°K; p = 0:7497+ 3-493. 10-* P—3-431.10-7 P?+ 
+0-01713.10-§ P® 
717-31°K; p = 0:8079+2:521.10-* P—2-031 . 1077 P?+- 
+0-007995 .10-® P® 
O;; 90:07°K: p= 1-1419+ 1-969. 10-4 P—1-233 P?. 1077+ 
+0-4682.10-1° P® 
77-35°K: p= 1-1997+ 1-608. 10-* P—0-8078 P?.10 7+ 
+0-298.10-1° P® 


DIE DICHTEN VON GESATTIGTEM, 
FLUSSIGEM WASSERSTOFF 


D._E. DILLER 
L. A. WEBER 


R. D. GOODWIN 

H. M. RODER 
Dre Dichte von Parawasserstoff von 17 bis 32°K wird durch Extra- 
polation der neu bestimmten Isothermen der komprimierten 
Fliissigkeit auf die entsprechenden Dampfdriicke erhalten. Die 
Ergebnisse werden mit friiheren Werten fiir den Bereich von 14 bis 
20°K kombiniert und mit normalem Wasserstoff mittels der Dichte- 
Temperatur—Beziehungen verglichen. 


DIE VERFLUSSIGUNG KLEINER 
WASSERSTOFFMENGEN 
B. YATES F. E. HOARE 
Es wird ein Verfahren beschrieben, das beim Entwurf eines Wasser- 
stoffverfitissigers fiir eine festgelegte Verfliissigungsmenge Anwend- 
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observées de la résistance sont beaucoup plus importantes lors de la 
transformation des alliages que lors de la transformation du lithium 
pur. 


UN INDICATEUR A NIVEAU D’HELIUM 
LIQUIDE 


A. E. ROVINSKII 


L’mwEE d’utiliser la transition du tantale pur état normal—état 
supraconducteur, a 4,3° K, pour mesurer le niveau de l’hélium 
liquide fut proposée par Feldmeier et mis en application pour la 
premiére fois par Rasor, dans un appareil pour la mesure continue 
de la quantité de liquide contenue dans un vase clos. Pour les 
indicateurs de niveau liquide a é/ément sensible en forme de sonde 
il est souhaitable de réduire les dimensions et la masse du tantale 
immergé, il est par suite, recommandé de changer le circuit 
électrique de l’instrument, de facon a conserver une sensibilité 
suffisante. 


ung fand. Es wird bei einem Bedarfsfall eingesetzt, bei dem 5 |. 
fliissigen Wasserstoffs ohne Verwendung eines Verdichters zu 
erzeugen waren. Konstruktive Einzelheiten des Reinigungs- und des 
Verfliissigungssystems einer zugdnglich und einfach gebauten 
Maschine sowie Leistungsziffern werden angegeben. 


TOLERANZEN BEIM BETRIEB DES 
CROWE-SPEICHERELEMENTES NACH 
DEM KOINZIDENZSTROMVERFAHREN 


R. V. PEACOCK 


Diz Arbeitsweise eines Crowe-Speicherelementes in Koinzidenz- 
stromschaltung wird beschrieben. Der Zusammenhang zwischen 
den Toleranzen in den AnordnungsgroBen und dem Betriebsstrom 
wird berechnet in Abhangigkeit vom Verlauf der adiabaten Wider- 
standskurve tiber dem Strom in der Querschiene (cross-bar). 


EIN VERFAHREN ZUR MESSUNG DER 
WARMEKAPAZITAT VON KONDENSIERTEN 
GASEN UBER DEREN SIEDEPUNKT 


A. V. VORONEL P. G. STRELKOV 


EIN Verfahren zur Messung der Warmekapazitat unter Wirbelung 
bei Driicken bis zu einigen Zehnereinheiten von Atmospharen 
wurde entwickelt. Die Warmekapazitat des Kalorimeters mit einem 
Inhalt von ca 105 cm* und einer Belastbarkeit bis zu 70 at lag 
zwischen 13-20% der gemessenen Menge. Der bei niedrigen Tem- 
peraturen und hohem Druck arbeitende Magnetwirbler fiihrt zu 
keiner bemerkenswerten Erwarmung und setzt die Zeit bis zur 
Herstellung des Gleichgewichts um das 1,5 bis 2-fache herab. 
Versuchsmessungen der Warmekapazitaét von Stickstoff im kriti- 
schen Gebiet wurden durchgefihrt. 


DIE MESSUNG DER WARMEDEHNUNG VON 
METALLEN BEI TIEFEN TEMPERATUREN 


K. ANDRES 


En Apparat zur Messung der Warmedehnung bei tiefen Tempera- 
turen wird beschrieben. Mit Hilfe eines optischen Verfahrens ist es 
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imOglich, Langenanderungen in der Gr6Be von 1 A (10-8 cm) bei 
Proben von 7 cm Lange zu bestimmen. Es werden Ergebnisse 
angegeben fiir den Gitteranteil und den elektronischen Anteil an 
der Warmedehnung von Al, Pb, Ta, Nb, Mg, und W. 


KOHLEWIDERSTANDSTHERMOMETER 
FUR TIEFE TEMPERATUREN 


N. N. MIKHAILOV A. YA. KAGANOVSKII 


Dre Technik der Herstellung von Kohlethermometern aus Oelpech 
zur Messung tiefer Temperaturen mit hoher Empfindlichkeit und 
befriedigender Reproduzierbarkeit wird beschrieben. Der Tempera- 
turverlauf des Widerstandes dieser Thermometer wird durch die 
Hrenntemperatur bestimmt. 


EIN WASSERSTOFFVERFLUSSIGER MIT 
EINEM AUSSTOSS VON 50 L/h 


A. G. ZEL’DOVICH YU. K. PILIPENKO 


/ARBEITSWEISE und Aufbau eines Verflissigers fiir 50 1./h Normal 
oder 30 1./h Para-Wasserstoff werden beschrieben. Es wird auf die 
Moglichkeit hingewiesen, die Anlage in einem Ka§4ltekreislauf zu 
verwenden. In einer Tabelle wird die Verfliissigungsleistung angege- 
ben, desgleichen die gemessenen Verluste. 


STANDANZEIGER FUR FLUSSIGES HELIUM 


A. E. ROVINSKII 


Die Idee, die Sprungtemperatur von reinem Tantal bei 4,3°K zur 
Bestimmung des Fltissigkeitsspiegels von Helium zu verwenden, 


Sommari 


LA VARIAZIONE DELLA DENSITA’ 
DELL’AZOTO E DELL’OSSIGENO LIQUIDI 
IN FUNZIONE DELLA PRESSIONE 


A. VAN ITTERBEEK O. VERBEKE 


SONO state eseguite misure sulla variazione del peso specifico 
dell’azoto e dell’ ossigeno liquidi fino a 850 Kg/cm7. I risultati possono 
essere rappresentati analiticamente dalle seguenti equazioni (p in 
g/cm’, P in Kg/cm?): 


N,: 90,26°K; p = 0,7497+ 3,493. 10-4. P—3,431. 10-7. P2?+ 
+0,01713. 10-8. P 


77,31°K; p = 0,8079 + 2,521 .10-*. P—2,031.10-7. P?+ 
+0,007995.10-®. P? 


O,: 90,07°K; p = 1,1419+1,969. 10-4. P— 1,233. P?. 10-7 + 
+0,4682. 10-9, P3 


77,35°K; p= 1,1997+ 1,608. 10-4. P—0,8078 10-7. P?+ 
=+- 05298) 10:20 Ps 
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stammt von Feldmaier und wurde erstmalig von Rasor in einem 
Apparat angewendet, mit dem die Heliumfiillung in einem geschlos- 
senen Behilter fortlaufend gemessen wurde. Bei einem Standanzeiger 
in Form einer Sonde ist es wiinschenswert, Abmessung und Masse des 
eingetauchten Tantalfiihlers zu verringern, weshalb man zweckma- 
Bigerweise einige Anderungen an der elektrischen Schaltung des 
Instruméntes vornimmt, um eine ausreichend hohe Empfindlichkeit 
zu gewahrleisten. 


DER EINFLUSS DER 
MARTENSITUMWANDLUNG AUF DEN 
ELEKTRISCHEN WIDERSTAND VON LITHIUM 
UND NIEDERPROZENTIGEN 
LITHIUM-MAGNESIUM-LEGIERUNGEN 


J. S. DUGDALE D. GUGAN 


Der Einfilu8 der Martensitumwandlung auf den elektrischen Wider- 
stand von Lithium und niederprozentigem Lithium-Magnesium- 
Legierungen wurde untersucht. Der Widerstand von reinem Lithium 
bei 80°K fallt bei der Umwandlung um etwa 0,5%, der von Legie- 
rungen betrachtlich starker. Versuche unter Anwendung von 
Warmebehandlungen und mechanischen Deformationen ergaben, 
daB der eigentliche Restwiderstand bei dichtester hexagonaler 
Packung und bei kubisch flachenzentriertem Gitter jeweils etwa 10 
bis 15% hoher ist als der bei kubisch raumzentriertem Gitter. Die 
Versuche ergaben groBe Abweichungen von der Mathiesen’schen 
Regel in allen drei Gefiigearten bei den Legierungen und zumindest 
in der kubisch raumzentrierten Anordnung von ziemlich reinem 
Lithium. Diese Abweichungen erklaren die wesentlich groBeren 
Widerstandsanderungen bei der Gefiigeumwandlung in Legierungen 
verglichen mit den bei reinem Lithium beobachteten. 


LA DENSITA’ DELL’>IDROGENO LIQUIDO 
SATURO 


R. D. GOODWIN 
H. M. RODER 


D. E. DILLER 
L. A. WEBER 


Sono state ottenute le densita del para-idrogeno tra 17°e 32°K 
estrapolando le isoterme, determinate nuovamente, del liquido 
compresso alle corrispondenti pressioni di vapore. I risultati soon 
stati combinati con dati precedenti per l’intervallo da 14° a 20° Ke 
confrontati con lo idrogeno normale per mezzo delle relazioni tra 
densita e temperatura. 


UN LIQUEFATTORE DI IDROGENO DI 
PICCOLE DIMENSIONI 


B. YATES F. E. HOARE 


Si da una descrizione del procedimento adottato nel progettare un 
liquefattore di idrogeno per produrre una quantita di liquido 
desiderata. La procedura é applicata al caso ove siano richiesti 5 It di 
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idrogeno liquido senza |’uso di un compressore. Sono inclusi i 
dettagli costruttivi dei sistemi di purificazione e liquefazione di una 
macchina accessibile nonché facilmente costruibile, insieme alle 
curve di prestazione. 


LE TOLLERANZE ASSOCIATE AL 
FUNZIONAMENTO DELLA PILA DI CROWE 
NEL MODO CORRENTE COINCIDENTE 


R. V. PEACOCK 


St descrive il modo corrente coincidente di funzionamento nella pila 
di Crowe. E’ calcolata la relazione tra le tolleranze sui parametri del 
dispositivo e le correnti in funzione della forma della curva adi- 
abatica della resistenza funzione della corrente della barra di collega- 
mento. 


LA MISURA DELLA ESPANSIVITA’ TERMICA 
DEI METALLI ALLE BASSE TEMPERATURE 


K. ANDRES 


E’DESCRITTA un apparecchiatura per la misura della espansivita 
termica a basse temperature. Usando di una tecnica ottica, é pos- 
sibile di rivelare variazioni di lunghezza di 1 A (10-8 cm) in un 
campione lungo 7 cm. I risultati sono dati in vista dei contributi 
reticolare ed elettronico alla espansivita termica nei metalli Al, Pb, 
Ta, Mg, e W. 


UN TERMOMETRO A RESISTENZA DI CARBONE 
PER BASSE TEMPERATURE 


N. N. MIKHAILOV A. YA. KAGANOVSKII 
Sono descritte la tecnica di preparazione dal bitume di petrolio e le 
proprieta di termometri a carbone, per misure di basse temperature, 
con alta sensibilita e soddisfacente riproducibilita. L’andamento 
della dipendenza della temperatura dalla resistenza di questi termo- 
metri é determinata dalla temperatura di cottura. 


Pe3iome 


M3MEHEHUE IWJIOTHOCTU ¢KUKOLO ASOTA UW 
HKUCJIOPOTA KAK ®YHKUMA JABJIEHWUA 


A. WAN ITTERBEEK O. VERBEKE 


Visamepunoch u3mMeHeHMe yeIbHOrO BeCa MAKOTO az0Ta M 
*#KMKOrO KucHOpofa go 850 xr/cm?. Pes3ynpraTEt MooxKHO 
TiIpeACTaBAATh aHaIMTMYeCKM IyTeM Ciefywulero ypaBHeHuA 
(po Br/cm?, PB xr/cm2): 
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UN LIQUEFATTORE DI IDROGENO CON 
UNA PRODUZIONE DI 50 It/ora DI 
IDROGENO LIQUIDO 


A. G. ZEL’DOVICH YU. K. PILIPENKO 


St descrivono le fondamentali caratteristiche di funzionamento e 
costruzione di un liquefattore che produce 50 It/ora di liquido 
normale, ovvero 30 It/ora di para-idrogeno. FE’ illustrata la possibilita 
di funzionamento del liquefattore in ciclo di refrigerazione. E’ data 
una tabella di prestazione del liquefattore e inoltre i risultati sulle 
misure di perdita. 


L’EFFETTO DELLA TRASFORMAZIONE 
MARTENSITICA SULLA RESISTENZA 
ELETTRICA DEL LITIO E DELLE LEGHE 
DILUITE LITIO-MANGANESE 


J. S. DUGDALE D. GUGAN 


E’sTaTo studiato l’effetto della trasformazione martensitica sulla 
resistivita elettrica del litio e delle leghe diluite litio-manganese. La 
resistivita del litio puro a 80° K decresce approssimativamente del 
0,5% durante la trasformazione e quella della lega decresce note- 
volmente di pil. Eseguendo variazioni cicliche di temperatura e 
deformazioni meccaniche si trova che le resistivita residue intrinseche 
delle fasi h.c.p. e f.c.c. del litio e delle leghe magnesio-litio sono circa 
i110e15 % rispettivamente, maggiori della fase b.c.c. Gli esperimenti 
mostrano che ci sono grandi discordanze dalla regola di Matthie- 
sens in tutte e tre le fasi delle leghe ed in almeno nella fase b.c.c. 
del litio abbastanza puro. Queste discordanze spiegano come mai le 
variazioni della resistenza, osservate quando le leghe si trasformano, 
sono molto pil grandi di quelle osservate quando si trasforma il 
litio. 


UN INDICATORE DI LIVELLO PER 
ELIO LIQUIDO 


A. E. ROVINSKI 


L’1DEA di usare il fenomeno di transizione nello stato superconduttore 
a 4,3°K del tantalio puro per misurare il livello dell’elio liquido fu 
proposto da Feldmeier, e per la prima volta applicato da Rasor in un 
apparecchio per la misura continua della quantita di liquido in un 
recipiente chiuso. In un indicatore di livello del liquido, con I’ele- 
mento sensibile in forma di sonda, é desiderabile di ridurre le dimen- 
sioni e la massa dell’elemento di tantalio immerso, cosicché é 
vantaggioso variare un poco il circuito elettrico dello strumento al 
fine di mantenere una sensibilita sufficientemente alta. 


Nz: 90:26°K; p = 0:7497+ 3-493. 10-* P—3-431.10-7 P?+ 
= O-017 13210 P* 


ie Wake = 10;8079452-521 Ose — 2:03) Onis 
+0:007995.10-8 P® 


Oz: 90:07°K; p = 1-1419+ 1-969. 10°>* P—1-233.1077 P*+ 
+0-4682.10-1° Ps 


TU So eis =a 1997 1-608. On. P—O0:807 8.1 Omura 
+ 0-298. 1071° Ps 
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Pesrome 
NJIOTHOCTM HACbINJIEHHOrO #KU KOO 
BOJJOPO]IA 


D. E. DILLER 
L. A. WEBER 


R. D. GOODWIN 
H. M. RODER 


Ilonyyawtca mmoTHocta napapoyopona or 17° go 32°K 
NYTeM SKCTPaiONANMM HOrO-OMpefeNeHHEIX M30TepM CHKaTOH 
FKMAKOCTH JO COOTBETCTBYWUIMX TaBNeHM mapa. Pesysb- 
TaTH COCAMHAIOTCA C PaHbINMMM aHHbIMM B Mpefenax 14° 
70 20° KH wu cpaBHMBaioTCA C HOPMAbHEIM BOJOPOOM WyTeM 
COOTHOWEHN IOTHOCTh — Temnepatypa. 


O}KMWIKAHUE BOJOPO]A B MAJIOM MACHIITABE 


B. YATES F. E. HOARE 


OnucaHMe MeTOAa pM KOHCTPYKUMM OKMHKaTeIA BOMOpona, 
KOTOPHIt MpoMsBo Ma Obl Tpe6oBaHHy!o MOWHOCTh. Meroy 
MCHONIb3yeTCA WIA cyyad, KOrga TpeOyeTcA BEIxOR DO I 
*AMAKOPO BOAOpOAa 6e3 ynoTpeOseHuA KOMUpeccopa. IIpu- 
COC€HMHEHbI KOHCTPYHTMBHbIe jeTaM CMCTeEMBL OYMCTKM UM 
ORKWHAHMA MAWMHbL, KOTOpaA Oba Obl WOcTyNHaA M MO]- 
XOJAMAA TIA KOHCTPyYKUMM, Tak *Ke HaK M TexXHMyeCKHe 
3HadeHiA. 


AONVYCK TPM SKCIJIYATAUMM AKKYMYJIATOP- 
HOLPO JJIEMEHTA KPOBE B BUTE C COBITA EHUEM 
TOKA 


R. V. PEACOCK 


OGcvewkyaeTCH BU OKCHIyaTayum Cc CcOoBMayeHweM TOKAa 
aKKYMVIATOPHOrO daemeHTa Hpose. BxtyncneHo oTHOMeHne 
Me*KY JONYCKAaMM TapaMeTpOB yCTaHOBKU WM IpMBOHEIMM 
TOKAMM B WYHKUMM MopMbI ayMadaTMyecKoh KpPMBOK Compo- 
TUBJICHMA B 3ABMCMMOCTM OT TOKA ApMa. 


METOTJMRA WU3MEPEHWUA TEDJIOEMHOCTM HKOH- 
AEHCMUPOBAHHBIX YTA30B BbIMNE UX TOUHM 
KMMEHUA 


A. V. VORONEL P. G. STRELKOV 


PaspaboTaHa MeTO]MKa M3MepeHMA TeMNOeMKOCTM C Tepe- 
M€IUMBaHMeM TIP jaBIIeHMAX JO HECKONbKMX JeECATKOB am. 
Tenmo0eMKOCTh Kaslopumetpa, BbijepsxuBalomero 70 am B 
o6beme ~ 105 cm, coctapnaAna oT 13 qo 20% mamepnemont 
BeIMYMHbI. MarnutHad Mellianka, padotawljad Up Hu3sKOK 
Temieparype M NOBBILMeHHOM JaBJICHMM, He CO3,aeT 3aMeTHOTO 
IpuToKa Telma MU cOKpallaeT BPEMUA YCTAHOBIeEHMA paBHO- 
Becua B 1,5+2 pasa. IIpouspeqentt npoOuEre musmepenna 
TeNMOCMHKOCTH a30Ta B KpUuTMYeCcKOH OOmacTH. 


M3MEPEHME TEPMMUUECKOLO PACIIMPEHMA 
METAJIJIOB IIPM HUSKUX TEMITEPATYPAX 


K. ANDRES 


OnncpiBaeTcA alllapar JIA M3MepeHMA TepMMYeCKOrO pacuM- 
peHMA pM HM3KMX TeMMepaTypax. C MOMOMIbIO OnTMYeCKuxX 
Me€TONOB MO*KHO OOHAPYHUTb B OOpaslax WNMHOM B 7 CM. 
usmenenna jimmnnt 1A (10-8cm). Tonasannt pesynpratet 1A 
pelleTHKM M BHA] bl 9EKTpOHOB [LIA TepMMYECKOTLO pacuiMpeHuA 
B Metasmax Al, Pb, Ta, Nb, Mg u W. 
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YIOJIbHbIE TEPMOMETPHI COMNPOTMBJIEHMA JIA 
; HU3KUX TEMMEPATYP 


N. N. MIKHAILOV A. YA. KAGANOVSKII 


OnucaHbl TEXHONOrMA M3rOTOBIeHMA M CBOMCTBAa YTOJIbHBIX 
TepMOMeTPOB M3 He*TAHOTO eka AA M38MepeHMA HUSKMX 
TemnepaTyp, OOmagqaiwoulMx BEICOKOM YyBCTBUTeIbHOCTbIO M 
YHROBIETBOPUTeEILHOM BOCIPOM3BONMMOCTbIO. XapakTep TeM- 
epaTypHOM 3aBMCMMOCTM CONPOTMBIIEHMA 9TUX TEPMOMETPOB 
onpenesiAeTcA TemnepaTypow oOmura. 


BOJOPOLHbIM OMMKUTEb MPOMSBOMTEJb- 
HOCTbIO 50 «/« HKMJ[KOTO BOJOPOTA 


A. G. ZEL’DOVICH YU. K. PILIPENKO 
OnmcaH mpHHouMm padoTh M KOHCTPpyKUMA OMMWKUTCIA 
npo“sBoquTenbuocTbw 50 2/4 HOpMasbHOrO usm 30 1/4 napa- 
Bomopoga. OTMedaeTcCH BO3MO*HOCTH padoOTH OPKMWKUTeIA 
B WMKIe peppuxepatopa. IIpupeyena Tadmmua moKasateren 
padorEl OKMKUTeEIA, a TAKKE PesyIbTAaTH 3amMepa MOTepb. 


BJIMAHUE MAPTEHCUTHOLO TEPEXOJIA HA JJIEK- 
TPAYECKOE COMPOTUBJIEHUE JIUTMA UM CILIABOB 
MATHMA C MAJION“ TIPAMECbIO JIMTUA 


J. S. DUGDALE D. GUGAN 


Uccaeqopanocb BIMAHMe MapTeHCMTHOTO Nepexofa Ha 9JIeK- 
TpMueCKOe CONPOTUBIEHMe IMTUA UM CHIAaBOB MarHuA C MaJOl 
mpuMecbio uTMuA. ConmpoTuBmeHwe 4uMcToro MTNA pH 
80° K mpm mepexoze nonukaeTcA IpMOmM3NTeNbHO Ha 0.0%, 
B TO BpeMA KAK CONPOTMBIeHMe CNIAaBOB NOHMKAeTCA Oomee 
3HauMTeIbHO. IlocpeCTBOM OMBITOB 10 TepMMYeCHOM WMKIN- 
3aluuu M MexaHMyYecKOl fedopmMayuu OL0 HalifqeHo, 4TO 
BHYTPeHHMe OCTATOYHBIE COMPOTMBIICHMA KPUCTAIIMYeCKUXx 
a3 IMTMA M CMIaBOB: reKCarOHaIbHOL C MOTHOH ymaKOBKOK 
M KyOMYeCKONM rpaHelleHTPMpOBaHHOM He MpeBLIMAaoT Gombe 
uem Ha 10 uw 15% coorsetcrBeHHO BHYTpeHHee OCTaTOUHOe 
cOompoTuBeHve a3sbr © KyOMYeCKOM OObeMHOMeHTPUpOBaH- 
HOM pemeTKOH. OKCHEPMMEeHTHL MOKasaIM, YTO y BCeX Tpex 
(bas CuuaBOB M y YMCTOrO IMTMA M0 KpaitHelt Mepe B dase 
c KyOu4eCKOl OOLEMHOLCHPMPOBAHHOL pelleTKON CyWIeECTBYWT 
OonbummMe OTKIOHEHHA OT MpaBuna Maruecena. OTumMn 
OTKIOHCHHMAMM OOBACHAeTCA TOT akT, YTO MpM Mepexope 
B clijaBax HadmofqawTcA ropa3sqo OombuIMe MsMeHeHMA 
CONPOTMBIEHMA, 4YEM TE, KOTOPble BOZHMKAIOT UpM Mepexose 
B 4MCTOM JIMTUM. 


MHIMKATOP YPOBHA KYUTKOrO PEJIMA 
A. E. ROVINSKII 


Vigen ucnomb30BaHvA ABIEHMA Nepexofa uMcToro TaHTama 
B CBepxmpoBoyAujee CcoctonuMe mpu 4,3° AH qua usmMepenua 
YPOBHA ?KMAKOrO reMA Mpeqtoxena Denbameiiepom [1] u 
npUMeHena BrepBble CepHHom [2] B ycrpoiictrBe AA HeMpepHIB- 
HOrO M38MepeHMA KOUMYECTBA +KUAKOCTM B 3aKPBITOM cocyse. 
B MHJMkaTOpe ypoOBHA *KMAKOCTM C WaTUMKOM B Bue ulyna 
HKCTATEMbHO BCEMCPHO YMCHbIUMTh pasmep M Maccy norpy- 
*KACMOTO TAHTAIOBOrO |JIeEMeHTa, MOITOMY JIA COXpaHeHMA 
HOCTaTOYHO BEICOKOM YYBCTBUTENBHOCTH Wemecoobpa3HoO He- 
CKOJIbBKO A3MeHMTh 9eKTpM4eCKy!O CxeMy MpMGopa. 
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The Variation of the Density of Liquid Nitrogen 
and Liquid Oxygen as a Function of Pressure 


A. Van Itterbeek and O. Verbeke 


Institut voor Lage Temperaturen en Technische Fysica, Louvain 


Received 8 July 1961 


COMPLETING our systematic investigation on the P, p, T 
relations of liquefied gases'? some measurements have 
been carried out up to 850 kg/cm? on liquid nitrogen and 
liquid oxygen. 

Only the variation of (¢p/0P),with pressure is measured, 
after integrating twice with respect to pressure and using 
as integrating constants values of the specific weights p 
and (dp/0P), of earlier measurements at zero pressure, we 
obtained the density as a function of pressure. 


(g/cm3) 


Our apparatus is practically the same as that used for 
our earlier measurements; only an adaptation to the high | 
pressure region is made.!>? Q 


The results for nitrogen are represented in Table | and 
graphically represented in Figure 1. 


Table 1. Specific Gravity of Nitrogen 
as a Function of Pressure 


T (K) P (kg/cm?) p (g/cm’) 


0) 500 1,000 
90-26 842:3 0:9030 Po (kg £m?) 
648-3 08782 
496-2 0-8599 Figure 1. Specific weight of nitrogen as a function of pressure 


375-2 0:8423 
279-9 0°8232 
199-1 0-8073 


131-0 ~- 0-7901 cs 
0 0:7497 oS 
77-31 574-2 0:9030 


431-0 08882 
319-2 0:8705 
2924 FI 0°8558 


150-3 0-8419 
82-0 0:8274 1.24 
0 0-8079 
Q 
These two isotherms are evaluated as third degree 
curves with the aid of the least-squares method: 
90-26° K;p = 0:7497+ 3-493. 10-4P— 3-431.10~-7P?+ 
’ —8p3 
as "9 200 400 600 800 1,000 
77:31°K3p = 0:8079 +2:521.10~-4*P—2-031.10~7P? + (ae (kg/cm?) 
+0-007995. 10~8P? Figure 2. Specific weight of oxygen as a function of pressure 
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Table 2. Specific Weight of Oxygen The experimental results obtained for oxygen are 


as a Function of Pressure represented in Table 2 and graphically represented in 
Figure 2. 
T CK) P (kg/cm?) p (g/cm’) Earlier measurements up to 150 kg/cm? suggested a 
linear dependence of specific weight with pressure. 
90-07 851-8 1:2490 However, more accurate measurements showed a real 
662-9 1:2322 dependence of (dp/0P); with pressure. 
5160 1-2170 Going up to the high pressure domain it appears very 
400-0 1-2039 
308-2 1-1918 clearly a dependence of (dp/0P), from pressure. 
230-0 1-1808 We express our thanks to Miss Deprit of the Centre for 
162-0 1:1712 Electronic Calculations directed by Professor G. Lemaitre 
98-0 1:1597 of the University of Louvain, who very kindly consented to 


calculate the analytic equation representing our measure- 


iy el roe ments. We also express our sincere thanks to the Belgian 
527-0 1:2663. Ministry of Education for the financial support of this 
406-0 1:2537 research which was given to our laboratory as a division 
oes cee of the National Centre for Research at Low Temperatures. 


232-0 1-2338 
158-2 12229 


The calculated curves are: 
90:07°K; p = 1:1419+4+ 1-969. 10-4P— 1-233P?.10-7+ 


ied ait 1. VAN I A dV. O. Physica 26, 931 (1960) 
reed : ’ v . VAN ITTERBEEK, A., and VERBEKE, O. Physica 26, 
77°35°K; p = 1:1997+ 1-608. 10-*P —0:8078P?. 10-7 + 2. VAN ITTERBEEK, A., and VERBEKE, O. Cryogenics 1, 77 (1960) 
+ 0-298. 10-1°P8 3. VAN ITTERBEEK, A., and VERBEKE, O. Physica. To be published 
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The Densities of Saturated Liquid Hydrogen 
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THis laboratory is engaged in a programme for the deter- 
mination of thermostatic, thermodynamic, and transport 
properties of para-hydrogen, emphasizing high densities. 
Reported here are new data for the density of saturated 
liquid para-hydrogen from boiling point to critical point. 
A review of hydrogen properties! includes densities 
of normal hydrogen from triple to critical points and of 
para-hydrogen up to the boiling point. More recent 
determinations of vapour pressures and critical constants 
have been made for both hydrogen modifications.** 
Values of constants used in our computations are given 
in Table 1. To avoid confusion with amagat density, 


Table 1. Constants Used for Computations 


para normal 
Te = 32-9844) 33-180!) 
Tr = 13-803 13-947) 
o,X 10? = 38-2029") 38-3024 
o.X 10° = 15:2672%) 14-9370) 


p, the symbol o is employed for density in gram-moles 
per cubic centimetre. In following sections two approxi- 
mate methods for representing the data are mentioned 
briefly, while two more-accurate methods are given in 
greater detail. 


Experimental 

The Reichsanstalt method of Holborn and colleagues,° 
exemplified by recent work at the Ohio State University,’ 
was employed to obtain the PVT isotherms.® 

Saturated liquid densities are obtained by extrapolating 
forthcoming isotherms for compressed liquid para-hydro- 
gen to the corresponding vapour pressures.? The 20°K 
equilibrium hydrogen (0-21 per cent ortho-hydrogen), 
herein called simply para-hydrogen, was produced by flow 
over iron oxide catalyst while filling the densitometer, and 
was confirmed by vapour pressure measurements.’ The 
isotherms are characterized by reproducibility and accur- 
acy for pressure and for temperature of about 1 part in 
10,000. Platinum resistance thermometry, through cali- 
bration by the NBS temperature physics section, rests upon 
the NBS low temperature scale of Hoge and Brickwedde’® 
as revised in 1955 by subtracting 0-01°C from the former 
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temperatures for fixed points and throughout the range. 
For density, precision and accuracy are 2 and 10 parts in 
10,000 respectively, the latter figure arising from the several 
volumetric determinations. By accuracy is meant here the 
authors’ estimate of the probable total cumulative effect of 
all the calibration, adjustment, and observational errors. 

To achieve the extrapolations, compressed liquid 
pressures are represented by polynomials in the density 
along isotherms, employing least-squares methods with a 
high speed digital computer. The tedious smoothing 
contemplated for the temperature dependence of the 
coefficients of these polynomials has not been applied here 
as the effect would not be significant. By means of the 
computer, an iteration procedure is applied for extraction 
of the density corresponding to an assigned (vapour) 
pressure. Accuracy of the extrapolations is not simply 
established. On some isotherms, data points occur 
virtually at the coexistence boundary, while on others they 
may be removed within limits ranging from 7-5 per cent of 
the saturation density at 32° to 2:0 per cent of the satura- 
tion density at 17° K. An important criterion therefore is 
the self-consistency of the set of densities obtained. 

Results for para-hydrogen from 17° to 32° are presented 
in the lower portion of column 3, Table 2. The upper 
portion contains previously determined low pressure 
densities, the smoothed results selected being within the 
original experimental deviations.* Normal hydrogen den- ~ 
sities are given for comparison in column 2 of Table 3. 
The upper group is from reference 2, while the lower group 
is from reference 11 as adjusted in reference 1. 


Representation of the Data 

Four expressions for representing the data have been 
examined. 

(1) Anapproximate relation for the density of saturated 
liquid para-hydrogen up to 31°K is 


(c—o,) = —0-001810+0-0092597(T,—T)'3 ... (1) 


The simplicity of this expression renders it useful for 
some applications. In the range 14° to 31° maximum de- 
viation in 0 is —0-15 per cent at 22°K. The deviation at 
32° is —0-51 per cent. This equation clearly does not 
extend to the critical temperature. 

(2) A well known relation is derived from the Verschaf- 
felt equation and law of the rectilinear diameter, '* 


2o/o, = 1—(1—2¢,/o,)0+(1 — 6)” eee) 
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where 6 = (T-T,/(7.-T,) 


Subscripts t and c refer to triple point and critical point, 
respectively, and n = 3, or n=1. Equation (2) may be 
rearranged to the linear form 


(cGy —2)" =m a+b(T,—T)' os 


where a and b are constants. 

The data plots of equation (2a) for para-hydrogen are 
markedly non-linear. Percentage deviations of o calculated 
by equation (2) for n= are presented in Figure 1, 
relative to smoothed data from Table 2, column 4. 

(3) A linear plot for experimental data is a valuable tool 
for its examination. There is a useful function for saturated 
liquid densities, valid nearly to the critical point,'? 


(T.-T)"/(a-—o,) = A+ BT wits(3) 


where m = ,4y. For para-hydrogen in the range 14° to 32° 
the function on the left-hand side of equation (3) varies 
only 7 per cent, while density varies 50 per cent of the mean 
value. E. H. Brown has shown at this laboratory that 
equation (3) is equivalent to 


ea) 


(1—0)"/Z = 1+00 ... Ba) 
Table 2. Densities of Saturated Liquid para-Hydrogen 

T(K) wes | a (exper.) | o(cale.) | 100 = 
13-803 | 0-0382029 0-0381998 0-01 
13-990 | 0-0381286 0-0381232 0-01 
14-000 '  0-0778 | 0:0381191 

14-990 0:0377102 0:0377030 0:02 
15-000 i 0-1327 0:0376987 

15-990 0:0372633 0:0372631 0-00 
16-000 | 0-2130 | 0-:0372586 

16-991 | 0-0367931 0:0368012 — 0-02 
17-992 , 0:0362989 0:0363159 —0:05 
18-993 | 0:0357845 0:0358047 — 0-06 
19-995 | | 0:0352509 | 0-0352642 — 0:04 

| 

17-000 | 0:3250 | 0-0368288 0:0367970 0-09 
18-000 0:4758 | 0:0362733 | 0:0363119 -—011 
19-000 0:6726 0:0357945 0:0358010 —0:02 
20-000 / 09228 0:0352753 0:0352615 0:04 
20-268 1-:0000 |! | 0:0351115 

21-000 1:2341 | 0:0346992 0:0346898 0-03 
22-000 | 1:6142 , 0-0341069 0:0340821 0-07 
23-000 | 2:0714 | 0:0334638 0-0334330 0-09 
24-000 2°6135 0-0327268 0:0327363 — 0-03 
25-000 | 3:2487 0-0319869 0:0319835 0-01 
26-000 3:9854 0:0311607 0-0311635 —0-01 
27-000 | 48315 | 0-0302443 0:0302610 — 0-06 
28-000 | 5-7950 | 0:0292449 0:0292534 — 0:03 
29-000 6:8890 | 0-0281325 0-0281060 0:09 
30-000 8-1210 0-0267471 0-0267588 — 0-04 
31-000 |  9-5000 0:0250788 0-0250921 — 0-05 
32-000 11-0470 0:0227923 0-0227821 0-04 
32-400 0-0214352 

32-700 0:0199780 

32-900 0-0182832 

32-984 12-7700 0:0152672 
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where 
es (o ~~ o.)|(% ra a.) 


and « is a small arbitrary constant. For several cryogenic 
fluids the function on the left-hand side of equation (3a) 
varies linearly with T up to about 6 = 0-95. It then must 
fall sharply to zero with eventual infinite slope at 0 = 1. 
Errors of several per cent in o may be expected for 6 > 0-95. 
Figure 2 presents data from Tables 2 and 3 plotted in the 
co-ordinates of equation (3). Filled circles represent new 
results of this report. Table 4 presents values for A and B 
of equation (3) determined by least squares, and the values 
for « of equation (3a) computed from the relation 


Pet TARE _.. 3b) 

(4) Data in the co-ordinates of equation (1) indicate a 
sigmoid plot, hence a polynomial in (T,—T)'/*. For tri- 
fluoromethane'* and for nitrogen,'* closely spaced data 


right up to the critical point have been represented accu- 
rately by 


Table 3. Densities of Saturated Liquid Normal Hydrogen 


T (K) o (exper.) a (calc.) 100 a 
13-947 0-0383024 0-0383038 — 0-00 
13-990 0:0382860 0-0382861 — 0-00 
14-000 0:0382819 

14-990 0:0378680 0-0378653 0-01 

15-000 0-0378609 

15-990 0:0374230 0-0374236 —0-00 
16-000 0-0374190 

16-990 0:0369530 0:0369594 —0-02 
17-000 0:0369546 

17-990 0:0364610 0:0364707 — 0:03 

18-000 0:0364656 

18-990 0-0359500 0-0359552 —0-01 

19-000 0:0359498 

19-990 0-0354200 0-0354101 0:03 

20-000 0:0354045 

20-380 0-0352100 -0:0351889 0-06 
21-000 0-0348263 

21-990 0:0342340 0-0342178 0:05 

22-000 0:0342114 

23-000 0:0335549 

23-990 0-0328390 0:0328580 — 0-06 
24-000 0:0328506 

25-000 0:0320908 

25-990 0:0312540 0:0312737 — 0:06 
26-000 0-0312650 

27-000 0-0303590 

27-990 0-0293600 0:0293629 —0-01 

28-000 0:0293522 

29-000 0:0282131 

29-990 0:0269260 0:0269034 0-08 

30-000 0:0268889 

31-000 0-0252776 

31-990 0:0231420 0-0231499 —0-03 

32-000 0-0231238 

33-000 0-0190252 

33-180 0:0149365 


CRYOGENICS - DECEMBER 1961 


(per cent) 


Deviation ——> 


| 
14 18 22 26 30 34 
Temperature ——> (°K) 


Figure 1. Relative deviations of equation (2). Results calculated by 
equation (2) minus smoothed data of Table 2, column 4, derived from 
equation (4) 


(0-0) = 3 af(Te-TY? 4) 


Table 5 presents values for the constants of equation (4), 
determined by least squares. The calculated results from 
equation (4), with corresponding deviations in per cent of 
a, are presented by Tables 2 and 3. The density of para- 
hydrogen is less than that of normal hydrogen at all tem- 
peratures. Figure 3 is an outline for the para-hydrogen 
results. 


Table 4. Constants for Equations (3) and (3a) 


para normal 
A = +0-148 434 +0-145 985 
B = —0-000450 198 — 0-000 447 314 
a = —0-060717 — 0-061 563 


o 
fe) 


(T-T)°Y(o-q) — 
R 


@ 
a 
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Figure 2. Density function. Open circles of the upper curve for para- 

hydrogen are smoothed data from reference 2; filled circles are new 

data of this report. Data of the lower curve for normal hydrogen up to 

20-390° K as smoothed from reference 2; above 21° K from reference 11 
as adjusted in reference I 
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Table 5. Constants for Equation (4) 


para normal 
a,x10® = +6-2675345 + 6:646 6259 
a,x10® = +1-4973511 + 1-026 5796 
a,x 10? = —0-18306903 + 0-085 4904 60 
a,x10® = —0-0206931 81 — 0-067 0342 34 


Mr. Thomas Strobridge contributed generously both 
advice and operable programmes for the digital computer. 
He brought to our attention reference 14. Continuing 
encouragement from R. J. Corruccini and R. B. Scott has 
made possible the programme on hydrogen properties. 
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Figure 3. Outline of density-temperature relation for saturated liquid 
para-hydrogen 
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THE refrigerants available in a number of laboratories 
which have entered the low temperature field in recent 
years are often limited to liquid helium and liquid nitrogen. 
A small supply of liquid hydrogen would provide a 
valuable extension of the facilities available in these labora- 
tories, as well as providing a more economical alternative 
to liquid helium, to workers whose requirements do not call 
for temperatures below 10°K. For these reasons the 
design of a small and accessible hydrogen liquefier, of 
simple construction, is presented here in a form which will 
enable other workers to design liquefiers with outputs to 
meet their particular requirements. 

The definitions of the symbols used in the equations are 
collected in the Appendix. 


Basic initial considerations 


From a knowledge of the variation of the enthalpy of 
hydrogen with temperature and pressure it is possible to 
calculate the rate of liquefaction at an expansion valve 
situated at the end of a countercurrent interchanger 
carrying gas under known conditions from the equation 


A= Hie + Hee me) 


Using data from the work of Woolley, Scott, and Brick- 
wedde,! and assuming a 100 per cent efficient interchanger, 
operating with a fore temperature of 70°K, calculation 
showed that an average of about 0-68 1. of liquid could be 
expected from a 165 ft? cylinder, the pressure of which fell 
from 120 to 30 atm during the liquefaction process. Allow- 
ing for the gas passed before the onset of liquefaction 
eight such cylinders are used in producing approximately 
51. of liquid hydrogen. 


The purification system 


A schematic flow diagram of the liquefier is shown in 
Figure 1. The gas flows first through the dryer, which 
consists of a copper tube of 23 in. 0.d., 1 s.w.g. wall thick- 
ness, and length 18 in., containing silica gel. It then flows 
through the ‘cleaner’. This consists of a similar, though 
shorter, tube. This contains packed nickel gauze and is 
immersed in liquid nitrogen, its purpose being to remove 
impurities with relatively high boiling points. From here 
the gas passes through two purifiers containing a total of 
about 300 g of granular activated charcoal. These take the 
form of helices, immersed in liquid nitrogen, each helix 
being made from a 9 ft length of 3 in. o.d. copper tubing of 
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13 s.w.g. wall thickness. Plugs of nickel gauze are included 
at the entries and exits of the dryer and purifiers, in order to 
prevent pieces of silica gel or charcoal from being dis- 
lodged, while a further nickel gauze filter is included at the 
position shown, as an additional precautionary measure. 
The silica gel and charcoal are both reactivated before a 
run by being pumped and baked. 


Waste. 


Cylinders 


nitrogen 
Liquid 
hydrogen — % 


I. First interchanger 
II. Second interchanger 


III. Third interchanger 
lV. Final interchanger 


Figure 1. Schematic flow diagram 


Interchanger dimensions 


(1) The final interchanger. Assuming a rate of flow of 
8-5 m3/hr (i.e. 300 ft?/hr), measured at room temperature 
before liquefaction commences, the corresponding mass 
rate of flow would be 


0:2076 g/sec i 


Assuming that the high pressure gas was cooled to 70°K 
before entering the final interchanger, for a 100 per cent 
efficient interchanger the mass of gas passing back through 
the low pressure side would increase from 0-1609 to 
0-1900 g/sec as the high pressure fell from 120 to 30 atm. 

Largely for simplicity of construction all interchangers 
have been constructed according to the tube-on-tube form, 
i.e. the high and low pressure tubes are soldered together 
along their lengths, and coiled into helices, as indicated in 
Figure 2. The relevant details of the interchangers are 
summarized in Table 1. 
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Table 1. Interchanger Dimensions 


High pressure Low pressure 
tubes tubes 
Inter- Length : 
changer | Outer| Wall | Outer| Wall (cm) ee 
diam. | thickness | diam.| thickness 
(in.) (s.w.g.) (in.) (s.w.g.) 
Ist $ 26 d De, 300 | Copper 
2nd os 20 B 26 200 | 70% Copper- 
30% Nickel 
3rd $ 26 100 | Copper 
Final x 26 Be 26 200 | 70% Copper- 
30% Nickel 


The thermal exchange coefficients are calculated in the 
following manner. 

(i) The convection coefficient per unit length, Ay. 
In accordance with the definition of this quantity, 
1 a : 


hy a 0.10C,9'( 5 dim 


7 (3) 


Substitution of the appropriate numerical data into this 
equation, corresponding to the situation where the high 
pressure tube is carrying gas at 120 atm pressure, gives 


hy, = 6-233 x 10~* cal/deg. cm sec 


for the high pressure tube. The corresponding value for 
each of the low pressure tubes is 


hy, = 0-984 x 10~? cal/deg. cm sec 
(ii) The overall coefficient of heat transfer per unit 
length, U;. 
Since there are two low pressure tubes in parallel, the 
appropriate form of the expression for U, in the present 


case 1S 
1 1 1 


eee + TF 

Uy hy 2h, 

neglecting the very small effect of the metal wall separating 
the gas streams. Therefore 

U,, = 1-496 x 10~* cal/deg. cm sec YG) 


The length of the interchanger is calculated in this 
manner. Defining the efficiency of the interchanger as 


_ 94-8 


it may be shown that this is related to the overall coefficient 
of heat transfer per unit length by the equation 


a ag Ga Sin te l—e.m, Guna Cx 
my Cy—m, C, l—-e 


..(5) 


€ 


0, L = ..(6) 


Assuming an efficiency of 95 per cent and substituting the 
value of m, corresponding to a high pressure of 120 atm, 


0-1609 g/sec together with results (2) and (4) in equation 
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(6) gives L = 123 cm. Repeating this calculation using data 
corresponding to the situation where the high pressure is 
30 atm yields an interchanger length of 144 cm, while a 
length of 245 cm would be required for an efficiency of 
98 per cent. In practice a compromise length of 200 cm was 
used. 

The pressure drops along the interchanger are calculated 
using the following method. The pressures dropped along 
the various tubes of the liquefier may be estimated from 
the equation 

Ap ___—'16RTh, ~=dM 
Di 10°72, Ds pdt 


Applying this equation to the most unfavourable case for 
the high pressure tube, when this is carrying gas at 30 atm 
pressure, shows that the pressure drop will be about 0-2 
atm, while the corresponding pressure drop in the low 
pressure tubes will be about 0-03 atm. These are both 
acceptably low. 

(2) The second interchanger. An attempt has been made 
to preserve the cooling obtained in passing through the 
purifiers by causing the high pressure tube to travel from 
the second purifier to the liquefier through an evacuated 
tube. The high pressure tube leaves this outer tube at K 
(Figure 2) and enters the second interchanger. A tempera- 
ture of 150°K was assumed for the temperature of entry 
of the high pressure gas to the second interchanger, while 
the low pressure return gas was assumed to enter at 70°K 
and to leave at about 150°K. Substitution of the appro- 
priate data into the equations given earlier yielded a value 
of 200 cm for the length of this interchanger, for an 
efficiency of 95 per cent. Subsequent observations of the 
low rate of boiling of the liquid nitrogen surrounding the 
second purifier have strongly suggested that the cooling 
derived in passing through the purifiers is greater than was 
assumed, indicating that a shorter second interchanger 
would suffice. 

The highest temperature of emergence of the high 
pressure gas has been calculated as 103°K, and occurs 
when the high pressure is 120 atm. 

(3) The third interchanger. Assuming this temperature 
of entry to the single-current interchanger, and aiming 
for a temperature of emergence of 70-1°K, the dimensions 
of a suitable interchanger may be calculated from the 
conventional equation 


—hL 
6.—4, = (8,—9,) exp ee hs (7) 
Assuming that 4 in. o.d. tubing of 26 s.w.g. wall was to be 
used it was shown, with the aid of equations (3) and (7), 
that L = 98 cm. 

(4) The first interchanger. In order to extract the cooling 
obtainable from the gas emerging from the second inter- 
changer at a calculated temperature of 146°K, a further 
interchanger was included in the system before the cleaner. 
The returning low pressure tubes were arranged to share 
the vacuum jacketing tube with the ingoing high pressure 
tube, so as to preserve the cooling obtainable from the low 
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pressure stream before this entered the first interchanger. 
Application of the appropriate thermodynamic data 
indicated that this interchanger should have a length of 
about 300 cm for an efficiency of 95 per cent using the 
tubing listed in Table 1. 

It may be added that the calculated values for Reynold’s 
number were adequate for turbulent flow in every tube of 
every interchanger. The calculated overall pressure drops 
suffered during passage through the interchangers were 
approximately 1 atm in the high pressure line and 0:2 atm 
in the low pressure line. 


The design 


The design of the liquefier is reproduced, approximately 
to scale, in Figure 2. The vacuum jacketing tube, mentioned 
earlier, is soldered into a bush in the top plate A, and the 
vacuum jacket is continued a short distance into the lique- 
fier in the form of a | in. 0.d. 24 s.w.g. wall stainless steel 


Oe ae 2a 3.8 4 
Seale (tn.) 


Positions of the bushes through which pass the following tubes: 


: Vacuum tube carrying high and low pressure tubes. 

: Float entry. c: Liquid nitrogen filling tube. 

: Vacuum/pressure gauge tube. e: Gas thermometer tube. 
Siphon tube. g: Expansion valve control rod tube. 

: High pressure inlet tube. j, k: Low pressure return tubes. 


s>>agD 


Figure 2. Scale diagram of the liquefier 
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tube, from which the gas tubes emerge at K. The purpose 
of the stainless steel termination is to prevent the O-ring 
seal of the plate A from becoming unduly cold. After 
passing through the second interchanger, D, the high 
pressure gas passes through the third interchanger, F, 
located at the bottom of the pumped nitrogen bath, E. 
This bath consists of a tube of stainless steel, 4 in. o.d., 
16 s.w.g. wall, and 14 in. long. The level indicator consists 
of an aluminium wire supported by a float of foam poly- 
styrene. This is covered by a glass tube and a seal is made 
by using a hat packing at J, where a protective slotted 
brass tube is also supported. The liquid nitrogen bath is 
pumped at C, and a vacuum gauge in the pumping line 
gives an indication of the temperature of the pumped 
liquid nitrogen. From this bath the high pressure gas flows 
through the final interchanger, H, and thence to the 
expansion valve within the liquid hydrogen receiver at L. 
The expansion valve control rod consists of a 4 in. dia- 
meter stainless steel rod, sealed by metal bellows at the top, 
and controlled as shown. The helical spring B prevents the 
valve from closing during evacuation. With this design 
only low pressure gas enters the space between the control 
rod and its surrounding german silver tube, o.d. 3 in. and 
0-1 mm wall. The low pressure return gas follows the path 
shown. 

The liquid hydrogen receiver consists of a vessel of 22 
s.W.g. spun copper, M, of capacity | litre. This is surroun- 
ded by a copper radiation shield attached to the base of 
the liquid nitrogen bath by screws. The liquid hydrogen is 
extracted via a 4 in. o.d. german silver tube of wall thick- 
ness 0-1 mm, which extends from the bottom of the 
receiver to an extraction valve outside the liquefier. The 
tube is vacuum jacketed from its passage through the 
plate G to the tip of the extraction valve. Within the 
liquefier the vacuum jacket takes the form of a 1 cm o.d., 
0-1 mm wall german silver tube, while this changes to 
standard 3 in. bore copper tubing outside the liquefier. 
Mitred joints, similar to those described by Stout, were 
used at the corners of the siphon tube and its vacuum 
jacket, asbestos spacers being used to prevent the inner 
and outer tubes from touching. The design of the extrac- 
tion valve was based upon that of Fairbank.? A helium 
gas thermometer bulb is soldered to the top of the liquid 
hydrogen receiver and is connected to anexternal Bourdon 
gauge by 1 mm bore, 0-1 mm wall german silver tubing. 
This thermometer serves to indicate the time at which 
liquid hydrogen starts to collect in the receiver. The outer 
vacuum jacket is pumped via a 3 in. diameter tube N which 
feeds the throat of a diffusion pump. This high vacuum is 
shared with that surrounding the incoming and outgoing 
tubes, and with that surrounding the siphon tube. Detach- 
able vacuum and water connections are used in the ancil- 
lary pipework, so that it is easily possible to lower the 
outer jacket and radiation shield for any servicing which 
may be found necessary in course of time. If, for any reason, 
access to the second or third interchangers should be 


required at any time these may be lifted out through the 
hole sealed by the plate A. 
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Performance 


In cooling the cleaner, the purifiers, and the nitrogen 
bath within the liquefier, 201. of liquid nitrogen are used. 
During a run 251. are consumed. Recovery of the liquid 
nitrogen from the cleaner and purifiers at the end of a run 
indicates that, in all, approximately 35 1. are consumed. 
The time between the commencement of flow and the 
onset of liquefaction is 15 min. The intervals between the 
onset of liquefaction and the times of extraction of litres of 
liquid are approximately 17, 20, 23, 27, and 50 min. 

If a compressor were to be used with the machine, to 
maintain the input pressure at 120 atm, and if the un- 
liquefied return gas was to be recirculated the overall 
consumption of gaseous hydrogen would fall appreciably, 
whilst the output should rise to a constant value of just 
over 3 1./hr. 

The total cost of constructional materials, gauges, and 
valves used in the liquefier and purification system was 
rather less than £200. The construction was undertaken 
by a technician without previous experience in this some- 
what specialized field, assisted part of the time by one of 
the authors (B.Y.) and the total time taken, as nearly as 
can be assessed, was some three months. 


The authors wish to thank two members of Bradford 
Institute of Technology: Dr. D. L. Smare for his support 
in obtaining the facilities which enabled the construction 
and installation to be undertaken, and Mr. D. Ives for 
his conscientious workmanship. 
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APPENDIX 


Definitions of symbols, in order of their appearance in 
the text: 


H, = enthalpy per unit mass of inflowing gas (cal/g) 
H, = enthalpy per unit mass of outflowing gas (cal/g) 
Hy, = enthalpy per unit mass of liquid hydrogen (cal/g) 
E = fraction of incoming gas which is liquefied 
C, = mean specific heat at constant pressure (cal/g. 
deg.) 
7 = mean viscosity (poise) 
D = internal diameter of tube carrying gas (cm) 
dm/dt = mass rate of flow (g/sec) 
6,, = temperature at which high pressure gas enters 
interchanger (° K) 
6; = temperature at which low pressure gas leaves 
interchanger (°K) 
6% = temperature at which low pressure gas enters 
interchanger (° K) 
L = length of interchanger (cm) 
my = mass rate of flow of high pressure gas (g/sec) 
my, = mass rate of flow of low pressure gas (g/sec) 
H = mean specific heat at constant pressure of high 
pressure gas (cal/g. deg.) 
C, = mean specific heat at constant pressure of low 
pressure gas (cal/g. deg.) 
Ap = pressure (atm) dropped along a tube of length 
AL(cm) 
R = gas constant (erg/mole. deg.) 
T = mean absolute temperature (°K) 
dM/dt = rate of flow (mole/sec) 
p = mean pressure in tube (dyn/cm?) 
6,, = temperature of liquid nitrogen bath (°K) 
6, = temperature of entry of gas to third inter- 
changer (°K) 
6, = emergent temperature of gas from third 
interchanger (°K) 
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THE Crowe cell and its mode of operation have been 
described by a number of authors.'~* This paper considers 
the tolerances on the device parameters and drive currents 
which are associated with the coincident current mode of 
operation. 


Mode of operation 


The two binary states of a Crowe cell are represented by 
the sense of the persistent supercurrents which may be 
made to flow round the D-shaped holes in the thin super- 
conducting film (Figure 1(a) and (b)). 


EN ge “~~ T 
/ ‘ 4 ‘ 
r Piers A 

ANE wi \e 


(a) (b) 


Figure 1. (a) ‘I’ stored; (b) ‘0’ stored 


Figure 2 illustrates a drive wire configuration for coin- 
cident current operation. Figure 3 shows graphically the 
relationship between the drive currents (solid lines) and 
the currents circulating in the film (dotted lines). Unity 
coupling is assumed between the drive wires and the 
crossbar of the ‘D’s, which may be shown theoretically to 
be a good approximation. 

Let us assume that there is a circulating current /, 
flowing in the cell representing a ‘1’ stored (Figure 3). 
Now it is necessary that the half read pulses shall have no 
destructive effect upon the stored information. Thus, if 
we apply a half read pulse Jp/2 along one of the drive 
wires of Figure 2, the sum of the current induced plus the 
circulating current must not be high enough for any 
appreciable resistance to be induced in the crossbar. Let 
the current at which the resistance becomes appreciable 
be Jc ;; then we may write the inequality 


I, +]p/2 < Tos eh) 


It is necessary, however, that when a full read current 
pulse is applied to the cell the induced current plus the 
stored current should be sufficient to induce resistance in 
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the crossbar which will dissipate energy and cause the 
circulating current to decay. The wattage dissipated in the 
now resistive bar by this current therefore decreases with 


Ccrt a 
X' wires 


’ . 
“Y’ wires 


Figure 2 


time. The heat produced is dissipated away at a rate 
dependent upon the temperature of the bar. Thus, the 
crossbar will cool and become superconducting again after 
a certain time. If the current decays much more rapidly 
than the heat is conducted away, then the current falls 
very nearly to zero before the crossbar becomes super- 
conducting, and therefore the current trapped in when the 
bar becomes superconducting, /;, is almost zero. How- 
ever, it may be shown that if the heat is conducted away 
quickly, then appreciable currents can be trapped in when 
the crossbar becomes superconducting. 

Now the resistance induced in the crossbar is a function 
of the drive current and, thus, the rate at which the 
circulating current decays is dependent upon the drive 
current. Thus, /;, one of the parameters of the Crowe cell, 
is a function of the current being used to switch it; this 
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Figure 3. Coincident current mode of operation 
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clearly complicates the operation of the cell. However, in 
this analysis we shall make the simplifying assumption 
that if the sum of the induced current and the stored 
current exceeds a certain value J, », and hence a sufficiently 
large resistance is induced in the crossbar, then the trapped 
current is fixed at a steady value J; and that any lower 
value of this sum results in incomplete switching, and an 
unpredictable value of /;. 

This assumption is clearly an oversimplification and, in 
fact, the device is likely to have tighter tolerances in 
practice than the foll6wing analysis predicts. 

Using the notation discussed above we may write the 
inequality which must be obeyed if a full read pulse J is 
to switch the cell 

Ti+lIp > Ie2 weer 2) 


When this read-out pulse is removed, sometime after the 
crossbar has become superconducting again, only rever- 
sible changes take place, subject to certain limitations, to 
leave a circulating current of Jp —/z, which represents a 
0’ stored. 

As no information may be written in by a half write 
pulse we may write a third inequality 


(RES Be Oe oe mes) 


and we may express the necessary condition that writing 
takes place when Jy is applied by 


TIp-—Ty+Iw > Ic2 ... (4) 


When this write pulse is removed we have a ‘1’ stored 
again. Thus 
Iw-Ty = I, ED) 


Now if we write Ic. as I¢,;+AlJ¢ and apply pense 
(5) to inequalities (1) to (3) we have 


Iwt+Tpl2 < Igitlt sek) 
Iptlw > Ioith+Ale C7) 
IptIw/2 < leith wy sD) 
Adding inequalities (6) and (8) 
IptIw < 4o1+44)/3 49) 


Thus we have a number of inequalities which have to be 
obeyed if the Crowe cell is to operate in a coincident 
current mode. We are therefore in a position to consider 
the tolerances on the drive currents, the device current 
parameters, and the relationship between them. 


Tolerances 

Let +100¢ per cent be the tolerance on the cell para- 
meters, and + 100g per cent be the drive current toler- 
ances. Now if the device is to operate at all in the coin- 
cident current mode, the following inequalities modified 
from inequalities (7) and (9) must still be valid 


(1—g)\UgtIw) > d+0de1+4+Ale). ... (10) 
(l+qg)Upt+lw) < 40 -9)Ue14+1n)/3 seal) 
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Thus, to find the limiting tolerances we must write these 
inequalities as equations, and dividing equation (10) by 
equation (11) we have 


(1—g)/(1+q) = 3114+9(14+1/n)/4U1—2d ... (12) 


where n= Ue +4q)/Ale 


Tolerance (per cent) on drive currents = 100g 


Tolerance (per cent) on device 
parameters Jc and J; = |00¢ 


Figure 4 


From equation (12) it is clear that if n equals 3 there is 
no tolerance at all on either Ic, 1, Iw, Jp, or Jy, and if nis less 
than three the cell cannot be used in a coincident current 
store. 

The relationship between the tolerances on the drive 
currents and the tolerances on the device parameters given 
by equation (12) is plotted in Figure 4. 


Comments on the resistance against current transition 


It is evident from Figure 4 that the value of has a most 
marked effect upon the tolerance question. Let us therefore 
try to get some picture of the type of resistance against 
current transition which will be needed for the device to 
operate at all in the coincident current mode. 

As the inductance of the cell for the original Crowe 
configuration is about 1 muH, andina working storea cell 
could conceivably spend as much as half of its time under 
the influence of half read pulses, then the resistance which 
is just tolerable for one hour of storage would be about 
10—'3 Q. Thus, Jc ; is the current at which R is approxi- 
mately 10~'3Q. As an order of magnitude Jc » may reason- 
ably be assumed to be the current through the crossbar at 
which R is of the order of 0-1 Q. Therefore AJ, is the 
current change required for the resistance to rise from 
about 10~!3 Q by a factor of more than 10'7, on the graph 
of R against J taken with short pulses. (This adiabatic 
curve is considerably broader than the isothermal d.c. 


graph.’) 
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The significance of this is made clearer when it is realized 
that there have been no adiabatic resistance against 
current characteristics published in the literature which 
have a transition width which would satisfy the inequality 
(Ic. 1+14)/Alc > 3, which is necessary for coincident 
current operation, even making the most favourable 
assumption that J; = I¢ 1. 


Conclusions 


The tolerances which are to be expected in operating 
a coincident current matrix of Crowe cells are extremely 
tight, even when the most favourable assumptions that J; 
can have a discrete value and that the coupling factor 
between the drive wire and the crossbar is the same for 
each cell in the matrix are made. In fact, unless the 
parameter ‘n’ of the adiabatic R against J curve is greater 
than three the device cannot be used in a coincident 
current storage matrix. 
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collaboration with Mr. N. D. Richards, and the author 
would like to express his thanks to Mr. Richards for his 
help and for many useful discussions concerning this toler- 
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and Dr. Lock, both of the Royal Radar Establishment, 
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Kostryukov et al.' developed a method of measuring the 
heat capacity of condensed gases at pressures not more 
than 1-2 atm. The present work extends this range of 
attainable pressures to several tens of atmospheres. This 
improvement in the method makes it possible to measure 
the heat capacity of the liquid phase above the boiling 
point, and also of two-phase systems up to the critical 
point. 

For elevated pressures it is very difficult to construct a 
calorimeter, which is light enough for its heat capacity, to 
be small compared with that of the specimen, while at the 
same time being sufficiently strong to stand the pressure 
without distortion. The construction shown in Figure 1 
has turned out to be very successful. The calorimeter is 
made of two tubes (1) and (2), pressed out with the domes, 
of stainless steel with wall thickness 0-18 mm. One tube is 
fitted over the other to a depth of ~ 5 mm and soldered. 
At the join a wide (20 mm) ring (3) is put on, made of the 
same tubing, and also soldered. A thin walled steel tube 
(4), of diameter 5 mm, passes along the vertical axis of the 
calorimeter to take the platinum resistance thermometer 
(6) wound on the brass sleeves (5) at the ends of the 
calorimeter and soldered. With this construction the 
calorimeter, sealed with soft solder and of internal 
volume ~ 105 cm?, stood up to above 70 atm without any 
noticeable deformation. The calorimeter weighed ~ 55 g 
and its own heat capacity was ~ 25 J/°C at room tempera- 
ture. For the measurements of the heat capacity of nitro- 
gen (see below) the heat capacity of the calorimeter 
amounted to 13—20 per cent of the total amount measured. 

The calorimeter heater (7) was a 0:1 mm diameter 
constantan wire threaded through a copper capillary 
(external diameter | mm) formed into a spiral along the 
inner surface of the calorimeter. The capillary containing 
the heater was filled with BF4 paste which was then 
polymerized. 

In many calorimetric experiments the fluid being investi- 
gated is stirred during the measurements. In particular, 
thermodynamic equilibrium of the system is very slowly 


+ Received by PTE Editor 28 October 1959: Pribory i Tekhnika 
Eksperimenta No. 6, p. 111 (1960). 
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established in measurements on the heat capacity of 
two-phase systems, and also near the critical temperature 
of liquids. Stirring considerably improves this and is even 
essential, in principle, for measurements in the critical 
region.” 


Figure I 


In the calorimeter described, stirring is achieved by the 
magnetic stirrer (8), consisting of two stainless steel disks 
fixed to the end of a hollow cylinder of permalloy strip. 
The stirrer is placed on the inner tube (4) as axis and is put 
into to-and-fro motion by switching a current through 
coils wound on the vacuum jacket at heights corresponding 
to the ends of the calorimeter. There are holes in the stirrer 
disks, which are displaced relative to one another, as 
shown in Figure 1, to improve the mixing. In order that 
the operating conditions of the calorimeter should not 
change under the action of the stirrer, the current in the 
coils is not switched off completely, but is only switched 
from one to the other. The method of thermal insulation 
and also the system for filling the calorimeter with gas 
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were in principle the same as described previously.! ? The 
development required only a high pressure valve.* 
The heat capacity of nitrogen in the critical region was 


(J7°C /mole) 


7—-> (%K) 


A Without stirring 


Figure 2. Temperature dependence of the heat capacity C, for nitrogen 
(at a density p = 0-3 g/cm?) 


© With stirring 


measured in order to try out and study the method devel- 
oped. The measurements were made on technical nitro- 
gen, containing ~ 1-5 per cent oxygen. The heat capacity 
was measured at constant volume, close to the critical 
value, and at various temperatures up to the critical. The 
results are shown in Figure 2. It can be seen that the 
temperature dependence, on crossing the curve of coexis- 
tence of two phases, has a jump similar to that found for 
ethylene, oxygen, and carbon dioxide.>’ We were not 
concerned here with a special study of this discontinuity. 
The heat capacity obtained for the liquid phase is close to 
that given in the literature. The form of the discontinuity 
found on crossing the two phase curve does not contradict 
theoretical ideas.® 
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Some heat capacity values were obtained without stir- 
ring in order to find how the stirrer worked. It can be seen 
from Figure 2 that these points show no systematic 
departure from the curve. Evidently the working of the 
stirrer produces no appreciable heat influx into the calori- 
meter during measurements. The time to reach equili- 
brium varied from 50 min far from the critical point to 1-5 
hr in the critical region itself. On stirring this time decreases 
to 15-30 min depending on the vigour of stirring. This 
problem needs further study right in the critical region. 

The error for a separate measurement of heat capacity 
does not exceed | per cent. The scatter of points on the 
curve, except for the point at 103-4°K is not outside this 
limit. 


The authors thank V. N. Kostryukov for discussion and 
interest in the work. 
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Ir is well known that the thermal expansion coefficient 
of solids tends to zero at low temperatures with the specific 
heat. This fact is expressed by the Griineisen relation 


BV [Kc SIG 


where « is the compressibility, c the specific heat per mole, 
V the molar volume, and yg the Griineisen ‘constant’. 
This is actually a mean value over the different y¢ 
corresponding to different vibration frequencies w of the 
solid. yg; may therefore exhibit changes with temperature 
which depend on the frequency spectrum and ultimately 
on the forces between the atoms or molecules of the solid. 

At very low temperatures one may expect the Griineisen 
constant of non-metals to become temperature inde- 
pendent, so that the thermal expansion coefficient B of a 
Debye solid should tend to zero as T°. In metals, however, 
both the specific heat and the expansion cofficient should 
show an additional linear term due to the electrons (see 
below, under Theoretical considerations). 

Whereas the order of magnitude of 8 at room tempera- 
ture is about 10° per degree, it decreases to about 107° per 
degree in the region where the lattice and the electron 
contributions are comparable, that is, below 10°K. For 
this reason, only measurements above 10°K have been 
available so far (see, for instance, Nix and MacNair,! 
Bijl and Pullan,? Gibbons,*? Swenson,’ Fraser and Hollis 
Hallett®). Very recently, White® ’ has described an appara- 
tus which is capable of detecting relative length changes of 
about 10~? by a condenser method; he has also given first 
results on the thermal expansion of copper, iron, beryl- 
lium, aluminium, palladium, and chromium at the 
temperature of liquid helium. The present author has 
developed an apparatus using an optical technique and 
having about the same sensitivity. Preliminary results on 
aluminium and lead have already been presented.®? In 
this paper a description of the method is given, including 
some further results obtained with it. 


Experimental method 


(1) Principle of measurement of length changes. The 
working principle of the apparatus is the transformation 
of length changes into changes of light intensity. The light 
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is passed through two closely spaced grids, one of which is 
held fixed, and the other of which is attached to the free 
end of the specimen whose expansion is to be measured. 
One grid is regular while the other has an irregularity in the 
middle as shown in Figure 1. Such a grid system has been 
described by Jones.!° !! In position (a) no light can pass the 
lower half of the grids, in position (b), where the second 
grid is displaced with respect to the first by one line 
width, the situation is reversed. In position (c) the two 
halves are equally transparent. Thus the output of two 
differentially connected photocells placed behind the two 
halves of the grids is zero. When one grid is fixed and the 
other one attached to the specimen, changes in length of 
the latter will result in changes of the output of the photo- 
cells. This is a linear function of the grid displacement for 
not too large an intensity of illumination. 

(2) Cryostat construction. For the sake of stability, the 
grids and the specimen have to be close together and it 
becomes necessary to place the grids within the helium 
bath. The problem of leading the light down into the 
cryostat has been solved by the use of Perspex light pipes 
as shown in Figure 2. The light, coming from a Philips 
high intensity lamp (type 13113 C/04 for 8 mm movie 


(c) 
Figure 1. Working principle of the grids 
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projectors), is focused into a light pipe after passing a 
cold light mirror and an interference filter (in order to 
remove the red and infra-red components). In the light 
pipe itself, it is totally reflected a few times on the polished 
surfaces until it falls on to the 45 degree surface at the 
bottom, where it is totally reflected towards the grids. 
Behind the grids, the light from each half is collected 
separately and led to the cryostat head again by two other 
light pipes at the end of which the photocells are mounted. 
The grids and the specimen are located in a brass frame- 
work which itself is firmly attached to the lower part of the 
three light pipes. The whole system is surrounded by a 
vacuum. The german silver tubes going to the top of the 
cryostat are connected to the light pipes at their lower 
ends. In this way, the apparatus is quite insensitive to 
shock. If the lower part of the cryostat moves with respect 


Cold light mirror 


Flexible joint 


Movable grid 
Calibration axle 


Calibration mechanism 


Brass framework} 


Vacuum can» | 4 


Figure 2. The cryestat 
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to the top (for instance, due to the boiling of liquid helium), 
then the grids and the light pipes move together so that in 
first order the illumination of the photocells is not 
changed. 

The grid connected to the specimen is suspended with 
flexible joints in the housing so that it can move up and 
down with its plane parallel to the other grid. This latter 
can be fixed by three adjusting screws so that the spacing 
between the two grids is about 0-05 mm. In a first execu- 
tion, the grids were drawn on tracing paper and photo- 
graphically reduced to the actual size of 10 x 11 mm witha 
line width of 0:15mm. More recently a second pair of grids 
was kindly made for us by the firm Dr. E. Rist AG, 
Prizisionsteilungen, Stafa, Switzerland. These are made of 
optically flat pieces of glass (25 x 25 x 2mm), one side of 
which is coated with a layer of chromium. On this side, 
the grid is scratched with a precision scaling machine. This 
pair of grids, which has a line width of 0-1 mm, has two 
advantages over the previous ones. Firstly, due to the 
greater precision, the optical efficiency is increased and, 
secondly, the light is reflected rather than absorbed by the 
lines so that less heat has to be carried away from the grids 
at low temperatures. This is important because the grids 
are in thermal contact with the housing and the absorbed 
heat has to be carried away by thermal conduction via 
the housing and a copper link to the helium bath, thus 
restricting the lowest temperature that can be obtained in 
the housing. If, for instance, the helium bath is at a 
temperature of 1-5°K, the temperature of the housing 
and the grids rises to about 2:5°K when the light is 
switched on. 

The lower end of the specimen is attached via a calibra- 
tion mechanism (see below) to the housing. Between the 
upper end of the specimen and the grid there is a system of 
lever arms with flexible joints as shown in Figure 2. This 
amplifies the length changes of the specimen by a factor of 
three. The specimen has a length of 7:5 cm and can be of a 
thickness between one and about three millimetres; it is 
held by two clamps. Between specimen and clamps some 
paper or hard rubber is inserted to isolate the specimen 
thermally from the housing. This insulation makes it 
possible to increase the temperature of the specimen above 
that of the housing. Heating the specimen from 1-6 to 
10°K raises the temperature of the housing from 2-5 to 
2:6°K, which does not affect the precision of the measure- 
ment. 

(3) Calibration. The calibration is done by observing 
the voltage produced in the differential photocells when 
the lower end of the specimen is moved by a known 
amount. Figure 3 shows the working principle of the cali- 
bration mechanism. When the axle A leading from the 
top of the cryostat down to the housing H is turned one 
turn, the toothed wheel W turns by one tooth. The axle 
of the wheel W consists of a left thread on one side and a 
right thread on the other, so that if the wheel W turns, the 
nut N moves. The movement of the nut N is transferred 
to the lower end of the specimen S by means of the lever 
triangle T which itself can turn around the fixed point F. 
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The soft springs shown in Figure 2 tend to pull the speci- 
men upwards and to push the nut N towards the wheel W 
thus minimizing the clearance between nut and thread. 

One turn of the axle A corresponds to a displacement of 
the specimen of 3-5 x 10~* cm with an error of about 3 per 


Figure 3. Calibration mechanism 


cent due to uncertainties in geometrical dimensions. 
Usually, a few calibration turns are made in the same 
direction. The corresponding output changes are traced 
by a pen recorder; they equal each other within about 
3 per cent, so that the total uncertainty of the calibration 
is not greater than 6 per cent. 

(4) Temperature measurement. Since it was intended to 
detect electronic contributions to thermal expansion in 
metals, that is to measure the expansion coefficient below 
10°K, carbon thermometers have been chosen for the 
temperature measurement. Two 100 Q Allan Bradley 
resistors are soldered to the specimen either at the two 
ends or in the middle and at one end. Temperatures are 
calculated from the formula 


ro log R 
~ (A+ Blog R)? 


due to Clement and Quinell.!* This should be valid in the 
temperature range from 1-5 to 20°K to within 0-02°. The 
constants A and B are determined in each experiment by 
calibration in the liquid helium range. It is believed that 
with the present procedure the accuracy of the tempera- 
ture measurements around 10°K is about 0-1°. In the 
samples of aluminium, lead, magnesium, and tungsten of 
thickness 3 mm, the largest temperature differences 
between the two ends are observed at 10°K and amount 
to 0:1°K (the heater is attached at one end). In the thin 
sample of tantalum of thickness 0-5 mm, this temperature 
difference amounted to as much as 0°3° at 10° K due to its 
bad thermal conductivity and to insufficient thermal 
insulation from the clamps. 

In all experiments the helium bath is kept at a tempera- 
ture of 1-5°K. The sample, thermally isolated from the 
housing, is connected by a thermal link to the helium bath. 
It can be heated electrically while the temperature of the 
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housing remains practically constant. With the light 
switched on, the sample has a temperature of about 1-6°K. 
Often, a constant drift is observed on the recorder 
corresponding to about 10 A/min (1 A= 10~8 cm), but 
this does not affect the accuracy of measurement, since 
the thermal relaxation time is about 2 sec. The heater 
is switched on during 20 sec and then switched off again. 
Figure 4 shows a typical recorder trace. 


OA { min 


Figure 4. Typical recorder trace 


The linear expansion coefficient is determined by graphi- 
cal differentiation of the observed thermal expansion 
curve. In all samples measured so far the volume expan- 
sion is simply three times the linear expansion since they 
all had cubic structures except magnesium which was 
polycrystalline. 


Results 


The results in so far are of a preliminary nature as all 
measurements have been made on one sample of each 
element only. The reproducibility of the results from dif- 
ferent samples of the same element has not yet been 
investigated, so that we know nothing of the dependence 
on purity. 

In Figures 5 and 6 the volume expansion coefficient B of - 
the elements investigated, divided by the temperature, is 
plotted against 77. It is seen that B can be fairly well 
represented by the formula 


B = aT+bT? 


This is indeed the temperature dependence one would 
expect from the theoretical considerations discussed 
below. The observed constants a and b are given in Table 
1. The errors quoted are those estimated from each experi- 
ment and do not, as mentioned above, include the possible 
scatter between different samples. The rapid drop of B in 
the normal state of niobium below 4° K may be due to the 
presence of spurious superconducting inclusions at the 
ends of the sample. The magnet coil around the sample 
(made of superconducting niobium wire) produced only 
3,100 oersteds, so that at the end the field was perhaps not 
high enough to destroy superconductivity. While heating 
from 1-6 to 4°K, these inclusions are ‘pushed out’, giving 
a negative contribution to the expansion because the 
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volume of niobium in the superconducting state is greater 
than in the normal state. 


Table 1. Thermal Expansion Data of Some Metals 


Element ax 10'°(deg-*) bx 10"\(deg~) 

Al iiee 8 66+ 06 
Mg 18+10 30 + 6 
Ta (normal) 33+ 2 Mf ae 3 

Ww 0+ 2 3 + 06 
Pb (normal) 30+ 60 420 +60 

Pb (superconducting) 0+30 360 + 60 
Nb (normal) 75 +30 45 +10 


For lead and niobium the thermal expansion both in the 
normal and superconducting states have been measured. 
In superconducting lead, B seems to follow a T? law at the 
lowest temperatures, which would be in agreement with 
the theoretical expectation on the basis of the two fluid 
model for superconductivity (see end of next section). In 
superconducting niobium, 8 does not quite follow this 
law and we do not yet have an explanation for its 
behaviour. 


Theoretical considerations 
Using one of the Maxwell relations 
(OV/0T), = (0S/0p)y 


one can calculate the thermal expansion coefficient from 


20 40 60 80 100 
Pe (°k2) 


© Al @ Mg (_] Ta (normal) —) wW 


A. Nb (superconducting) 


Figure 5. The thermal expansion coefficient B, divided by T, for some 
metals plotted against T* 


© Nb (normal) 
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the pressure dependence of the entropy at constant 
temperature. The entropy of a metal is the sum of a lattice 
and an electronic contribution: 


S= §,+8, 


and the temperature dependence of the expansion co- 
efficient must also consist of two contributions due to the 


(10° Ydeg2) 


20 
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Figure 6. The thermal expansion coefficient B, divided by T, for lead 
plotted against T* 


lattice and the electrons. At very low temperatures c,, the 
lattice specific heat, varies as T° and in this region S, also 
varies as T°. This means that at very low temperatures the 
lattice contribution to thermal expansion f, varies as T°. 
At the higher temperatures S; is usually a more compli- 
cated function of temperature. There, the Griineisen 
parameter yg usually changes and various calculations 
exist to predict its temperature dependence (see, for 
instance, Barron,'? Blackman,'* Bijl and Pullan’). Cal- 
culations of the electron contribution have been given by 
Visvanathan!> and Varley.'® The entropy of the electrons 
equals y7, the electronic specific heat per mole. Therefore 
B., the electron contribution to the thermal expansion, 
also varies as T. In terms of y7, we calculate 


Be = K(y/V) eT 


where g stands for dlogy/dlog V and (y/V)T is the elec- 
tronic specific heat per unit volume. Since y is propor- 
tional to n(€), the electron density of states at the Fermi 
surface (see, for instance, Mott and Jones!’), g is also 
equal to dlogn(é)/dlog V. The linear term in the thermal 
expansion coefficient of metals is thus a measure of the 
relative change of the density of electron states at the 
Fermi surface with volume. In Table 2, the observed g 
values of the measured metals are tabulated and compared 
with White’s’ results. It should be noted that the room 
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temperature compressibilities have been used for the com- 
putation of our values of g. 


Table 2. The Quantity g= @logn(é)/dlogV 
for Some Metals 


g kx JO"? g 
TANG (this work) (cm?/dyn) (from White) 

Al Leff 0992 1-4 ESeEO? 
Mg 0-6+0:3 3-0 

Ta 1-32-02 0-5 1-4* 

W 0 +02 0:3 

Pb U7 fa 2 eas) BS} ila7f- 

Nb 211 0:5 


* Private communication. 


One can think of another more direct way of calculating 
the electron contribution to thermal expansion. The free 
electrons in a metal form a gas of particles which exerts a 
pressure on the metal walls. This internal pressure P can 
be calculated in a way similar to that in the kinetic theory 
of gases, if it is supposed that the Fermi surface is known. 
Its slight temperature dependence causes the electronic 
contribution to thermal expansion. For the special case 
of a Fermi surface with cubic symmetry, one can write 


Be i K(OP/0T)y 


where « is again the compressibility of the metal. A calcu- 
lation of (0P/0T)y in this case gives, under the assumption 
that the electron density of states at the Fermi surface does 
not vary too rapidly with energy 


Be «(OP/0T)y 


a 
KTV) —z k* N@logn(E)/9E)~ ¢] 


Here, N is the number of free electrons per unit 
volume, n(E) the electron density of states, and € the 
Fermi energy. Thus, the electron contribution to thermal 
expansion is also a measure of the first derivative of the 
density of states function at the Fermi energy if the elec- 
tronic specific heat yT is known. A detailed derivation of 
the above formula will be published elsewhere. 

If the metal becomes superconducting, the electrons 
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condense gradually into a state of zero entropy as the 
temperature decreases. In the picture of the two fluid 
model of Gorter and Casimir'*® the entropy of the electrons 
in a superconductor is written as 


S.= yT?/Te 


where yT is again the electronic specific heat per mole, and 
T, is the critical temperature of the superconductor. 
Inserting this in Maxwell’s relation, one obtains 


Pes ae K(y/V)\(T?/T2)(g — 25) 


where s is written for dlog7,/dlog V and the suffix s in 
B.s tefers to the superconducting state. For most of the 
superconductors, g—2s is negative, as Olsen and Rohrer”! 
have shown. The measurements on lead are in qualitative 


agreement with this statement. 


My thanks are due to Professor P. Grassmann, the 
Director of this Institute, for his constant interest and 
encouragement in this work and to Professor J. L. Olsen 
and Dr. H. Rohrer for much valuable advice and many 
discussions. [ am grateful to Dr. G. K. White for inform- 
ing me of his results prior to publication. This work was 
supported financially by the ‘Arbeitsbeschaffungsfonds 
des Bundes’. 
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TEMPERATURE measurement below 20°K presents well 
known difficulties mainly connected with the fact that 
thermocouples and resistance thermometers of metallic 
conductors become very insensitive in this region. One of 
the convenient methods of overcoming this difficulty is 
the use of carbon resistance thermometers which, with 
semiconductor conduction characteristics, have a negative 
temperature coefficient of resistance and increased sensi- 
tivity with decreasing temperature. During the last twenty 
years different variants of such thermometers have been 
published on numerous occasions. 

The simplest method is the application of a thin film of 
graphite or carbon black from the corresponding suspen- 
sions in water on a strip of insulating material, e.g. on 
cigarette paper.! Such thermometers are easy to prepare, 
but they are not sufficiently stable. Zavaritskii and 
Shal’nikov” have described a more reliable technique for 
preparing carbon film thermometers. 

Several types of radio resistor have found wide applica- 
tion at low temperatures in recent years,*’ and their pro- 
perties turned out by chance to be suitable for thermo- 
meters. 

Sharvin, who has prepared carbon thermometers from 
anthracite,® noted that the region of their greatest sensi- 
tivity can be moved in any direction by changing the 
temperature of firing. It was natural to suppose that this 
peculiarity is not a property of anthracite alone, but of a 
number of other organic substances with high carbon 
content. In this connection we have undertaken experi- 
ments to develop the technology of preparing sturdy 
carbon thermometers, starting from materials which 
allow greater reproducibility and uniformity of the 
thermometers than can be obtained by using anthracite. 


The technique of preparing thermometers 


The initial material is petroleum pitch, part of which is 
used without any treatment (in the form of brickettes 
bonded under pressure) and part is transformed into 
petroleum coke. Coking of the pitch is achieved by heating 


+ Received by PTE Editor 7 June 1960: Pribory i Tekhnika 
Eksperimenta No. 3, p. 194 (1961) 
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it in the absence of air at 700°C for 3 hr. The coke and the 
pitch are separately ground in a mortar and passed through 
a 625 holes/cm? sieve. It was found that to obtain the best 
results the powders had to be thoroughly mixed in the 
proportions of 23 per cent by weight of pitch and 77 per 
cent coke. The charge is then weighed out into portions 
which are pressed at room temperature in a hydraulic 
press at a pressure of 12,000 kg/cm”. For better bonding, 
each loading is kept under pressure for 5-10 min. The 
brickettes, obtained as rectangular parallelepipeds of 
dimensions 7x3x1 mm, are fired in a quartz tube 
filled with finely divided charcoal. The temperature is , 
raised over a period of ~ 1-5 hrand held at the given maxi- 
mum temperature for 1 hr. After this time the oven is 
turned off and the brickettes cool with the oven. 

Thin layers of copper are then deposited on to the ends 
of the brickettes by electrolysis from an acid bath and 
form cap-shaped contacts. Two 0-2 mm diameter copper 
wires are soldered at each end to the copper plated surfaces 
with a small quantity of tin, using a microsoldering iron. 
The thermometers are washed in spirits and covered with 
a layer of BF-2 adhesive which is polymerized in 2 hr at 
180°C. 


Results of the study of the thermometers 


The resistance of the thermometers was measured 
potentiometrically using a PMS-48 potentiometer and a 
M-25 galvanometer. After some preliminary experiments 
it was established that the optimal value of the measuring 
current was 25 mA. 

Preliminary experiments, carried out to choose the 
firing conditions, showed clearly (see Table 1) that raising 
the firing temperature leads to a reduction in the resistance 
of the thermometers. 

This relation can be seen at any temperature of measure- 
ment, but is most strongly marked at low temperatures. 
As can be seen from the table, there is considerable 
scatter in the resistance values within each group of 
thermometers which have been fired at the same tempera- 
ture. It is probable that this is due to the conditions of 
pressing (the brickettes not being exactly the same), and 
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mainly to small differences of firing temperature produced 
by non-uniformity of temperature inside the furnace. 


Table 1 
Thermo- | Firing Thermometer resistance (Q) 
meter | tempera- a 
number pee 
(CC) | 273°K | 77°K | 20:3°K | 4:2°K 

5-1 750 2-4 10-2 

5-2 750 2-7 12:6 

5-3 750 2:35 oe 

5-4 750 21 93 

5-5 750 2°5 12:2 

5-6 750 2-6 q°7: 

7-1 770 1-5 5°75 

7-2 770 1:3 fa 

13 770 1-65 6:2 Above 

7-4 770 1:75 EG 5,000 

7-5 770 1-87 68 

7-6 770 1:8 6-65 

8-1 780 1:7 65 

8-2 780 1:6 5-4 

8-3 780 2:0 8:1 

8-4 780 1-6 4-5 

8-5 780 1:62 4-8 

8-6 780 1:7 65 

9-2 790 1:26 S8/ 32:0 3,855 

9-3 790 1:17 2-7 71 959 

9-4 790 1:23 3-6 34+] 4,768 
80-1 800 0-9 1:77 6:7 145-8 
80-2: 800 0-7 1:37 41 47-7 
80-3 800 0-6 1:0 25 ie) 
80-5 800 0:45 0:66 1:3 66 
80-6 800 1-1 2-0 8-0 200 
81-1 810 0:55 0-95 25-7 
81-2 810 0:5 0:8 13-15 
81-3 810 0-54 1:0 33-5 
81-4 810 0:54 0-9 12°1 
81-5 810 0:5 0-8 8-75 
81-6 810 0-6 0-85 12:05 


Typical temperature variations of resistance for thermo- 
meters fired at three temperatures, 790°, 800°, and 810°C, 
are shown in Figure 1 (logarithm of the resistance as a 
function of the reciprocal of the absolute temperature). 

It is generally assumed that the temperature dependence 
of resistance of semiconductors in the intrinsic region can 
be described approximately by the equation 


p = cexp(W,/2kT) ert) 


where W, is the activation energy for the creation of 
carriers. 

As can be seen from Figure 1, this simple formula only 
gives the general character of the dependence for carbon 
thermometers, and is insufficiently exact for practical use 
in temperature measurement. Calculations show that 
between 20° and 2°K the interpolation formula proposed 
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by Clement and Quinnel? gives satisfactory results for our 
carbon thermometers: 


log R+K/logR = A+B/T ete) 


where K, A, and B are constants. 


a og 


0:1 0:2) 5 O36 0:4 0:5 0:6 


— (1/T deg. K") 
Figure 1. The temperature dependence of the resistance of thermo- 
meters which have had different thermal treatments. Thermometer 
80-2 was studied before and after 100 repeated cycles of cooling and 


warming from room to nitrogen temperature: 1, before the 100 cool- 
ings; 2, after the 100 coolings 


If these constants are determined from the resistance 
values at 2°, 4-2°, and 20-4°K, temperatures from 2° to 
4:2°K can be calculated with an accuracy of a few 
hundredths of a degree. To achieve greater accuracy a 
detailed calibration of each thermometer is essential, and 
the results can be tabulated or expressed as a graph 
T=f(R). Vroomen® considers it more convenient in 
practice to construct for each thermometer a graph of the 
corrections (‘deviations’) to function (2) which can be 
calculated with great accuracy and tabulated for any one 
thermometer of a given type. 

As can be seen from Figure 1, the curve log R = f(T~') 
becomes flatter the higher the firing temperature and 
corresponds to a shift of the region of a rapid increase of 
resistance in the direction of ever lower temperatures. The 
constant B in equation (2) which characterizes the slope 
of the curve and can be considered proportional to the 
activation energy rapidly decreases, as is seen from the 
curve in Figure 2 (constructed on the basis of studies of 
three thermometers Nos. 9-2, 80-2 and 81-4). 

It was important to examine the reproducibility of the 
thermometer calibrations. One must distinguish two 
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cases: (1) the reproducibility during one helium experi- 
ment and (2) the reproducibility after one or several 
cycles of cooling and heating up to room temperature. 
Experiments showed that during one helium experiment 
there is complete reproducibility, i.e. the errors in measur- 
ing the temperature are determined only by the accuracy 
of reading the helium vapour pressure if, as in our case, an 


790 800 810 820 


firing ar ee (CC) 


Figure 2. The variation of the constant B of the interpolation 
formula (2) as a function of the firing temperature of the carbon 
thermometers 


ordinary mercury manometer is used. On the other hand, 
differences of resistance after repeated cycles of cooling 
and heating are easily observed. The temperature varia- 
tion of resistance of thermometer 80-2 is shown in Figure 
1 before and after a hundred coolings and heatings from 
room temperature to nitrogen temperature. It can be seen 
from Figure | that the change in calibration after such 
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treatmént is small. In the worst case the ambiguity in 
measuring temperature in the range from 2° to 4-2°K did 
not exceed ~ 0-01°K. 


Conclusions 


A technique has been developed for preparing carbon 
resistance thermometers for low temperatures with suf- 
ficiently reproducible properties, using petroleum pitch 
as the only initial material. It is shown that the form of the 
temperature dependence of the resistance of the carbon 
thermometers is determined by the temperature of firing. 
Thermometers fired at 810°C are most convenient for 
measurements in the helium region. Thermometers fired 
at 790°C and 800°C are suitable for measurements in the 
hydrogen temperature region and in the intermediate 
range between hydrogen and helium temperatures. 
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THE use of liquid hydrogen bubble chambers in physical 
experiments has made it necessary to provide hydrogen 
liquefaction plants with an output of 50-200 1./hr. The 
hydrogen liquefier described was designed as an inde- 
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——-— Low pressure hydrogen 
—-:— Pumped nitrogen 


—— High pressure hydrogen 
+--+ Gaseous nitrogen 


1. Warm region heat exchanger; 2, Bath of liquid nitrogen, boiling 
at atmospheric pressure; 3, Heat exchanger between the nitrogen 
baths; 4, Bath of liquid nitrogen, boiling under reduced pressure; 
5, Cold region heat exchanger; 6, Throttling valve; 7, Additional 
collector; 8, Helium valve; 9, Reaction tube ; 10, Condenser coil; 
11 Draw-off valve for normal hydrogen; 12, Draw-off valve for 
para-hydrogen; 13, Draw-off valve for normal hydrogen from 
the additional collector; 14; Insulating Dewar 


Figure 1. Arrangement of liquefier 
eer re re ee ee ee 
+ Received by PTE Editor 19 April 1960: Pribory i Tekhnika 
Eksperimenta No. 2, p.185 (1961) 
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pendent plant or as the hydrogen section of a combined 
hydrogen-helium liquefier. The construction of the lique- 
fier provides for obtaining liquid para-hydrogen and for 
operating as a refrigerator. 

The arrangement of the liquefier is shown in Figure 1. 
Liquefaction is produced by a throttling cycle with liquid 
nitrogen precooling. The preliminary precooling takes 
place in two stages: by nitrogen boiling at atmospheric 
pressure, and by nitrogen boiling under vacuum. The 
special feature of the arrangement is the presence of an 
additional collector, 7, with the valve 8. The through cross- 
section of the valve is regulated by changing the pressure 
of helium in a bellows volume. 

For liquefying normal hydrogen the helium valve 8 is 
fully open and the liquid hydrogen runs over from the 
collector 7 into the lower part of the insulated Dewar 14, 
used as the collector of normal hydrogen. To obtain para- 
hydrogen, the draw-off valve 12 is opened. The liquid 
hydrogen passes from the additional collector into the 
reaction tube 9, containing a catalyst, and then into the 
coil condenser 10. The ortho-para heat of conversion given 
off in the reaction tube is partly taken off through its heat 
transfer surface. The remaining part of the heat produces 
evaporation of some of the converted hydrogen. Conden- 
sation and undercooling of the para-hydrogen obtained 
takes place in the coil condenser. The heat transfer is 
achieved by the normal hydrogen boiling in the Dewar. 
The quantity of para-hydrogen drawn off is regulated in 
such a way that a fraction of the liquefied normal hydrogen 
flows from the additional collector 7 into the Dewar 14, 
maintaining a fixed level there. The necessary temperature 
difference in the condenser coil is produced by increasing 
the pressure in the additional collector to 2-5-3 atm by 
using the helium valve. The reaction tube of the liquefier is 
made in the form of a coil out of 38 mm diameter copper 
tubing. The catalyst is ferric hydroxide, prepared in our 
laboratory on the basis of Buyanov’s studies.' The volume 
of catalyst was reckoned according to the accepted speci- 
fications? and amounted to 1,300 cm’. 

The inclusion of the additional collector with the helium 
valve in the liquefier construction makes it possible to 
recover the saturated vapour pressure in the liquefier 
after using the necessary latent heat of evaporation of 
hydrogen. The refrigerating capacity is doubled by using 
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the cold from the evaporating hydrogen (in the refrigera- 
tion cycle). In this case the liquid hydrogen is drawn off 
from the additional collector through valve 13 and the 
evaporated hydrogen returns to the lower part of the 
Dewar through the draw-off valve 11 or through a special 
siphon. It is convenient to use such a system of operation 
of the liquefier in order to maintain a constant thermal 
regime in liquid hydrogen bubble chambers. 

The liquefier is thermally insulated by a metal Dewar, 
filled with Mipor pumped to fore-vacuum pressure; the 
Mipor layer is 80 mm thick. Free space is left between the 
Mipor and the outer wall of the Dewar, to improve the 
pumping.’ There is no liquid nitrogen shielding. The 
liquefier apparatus is in hydrogen vapours within the 
Dewar, like a number of other plants.* The lid is at room 
temperature and the base at liquid hydrogen temperature. 
The absence of heat exchange by convection between 
parts of the liquefier at different temperatures is achieved 
by (1) the warmer parts of the apparatus being above the 
colder; (2) the temperature of all parts in any horizontal 
section being the same; (3) insulation of any tube for 
which the temperature does not correspond to that of a 
given section. 

The type of insulation used has several advantages: 
(1) leaks which frequently spoil the working of apparatus 
using the usual vacuum insulation are not dangerous in 
this case; (2) it is easy to achieve vacuum tightness of the 
Dewar because of the small number of joins; (3) it 
facilitates access to the heat exchange apparatus. 

The heat exchangers are made of bundles of soldered 
tubes. Since the heat flow along the tubes is very small, all 
the heat exchangers (including the heat exchanger of the 
‘cold’ region) are made of copper. The inner shell of the 
Dewar is made of stainless steel, the outer of ordinary 
carbon steel. The tubes for drawing off the liquid hydrogen 
and also the nitrogen siphons, have their own vacuum 
insulation. The usual level devices for low boiling liquids 
are used in the liquefier. These only give correct readings 
if the tube going to the lower part of the vessel is full of 
vapour, produced by boiling off the liquid. In our case the 
necessary heat supply is provided by the heat conductivity 
of the tubes and by the extra layers, which are insulated to 
prevent heat leakage. 

The performance of the liquefier is givenin Table 1. The 
experimental measurements showed that the refrigerating 
capacity of the liquefier (No. 2 in the table) is used in the 
following way: 


Q (kJ/hr) (per cent) 
On liquefying hydrogen 191 82 
As a result of non-regeneration 40-6 L725 
and pressure losses of the 
reverse flow 
On heat flow from the surround- Ie 0-5 
ings and the spontaneous 
ortho-para conversion 


Total 232:9 100 
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The heat flow from the surroundings is determined from 
the amount of hydrogen boiled off from the Dewar after 
the end of the operation. The loss resulting from non- 
regeneration is determined from the difference. If the 
heat flow through the insulation is referred to the parts of 
the Dewar inside which the hydrogen region apparatus is 
contained, i.e. to the region lower than the bath of pumped 
nitrogen, then it amounts to 1-4 kJ /m*. The consumption 
of liquid nitrogen is 1-3-1-35 1. per litre of liquid hydrogen. 
It takes 1-5 hr to cool the apparatus to the start of lique- 
faction, and this requires 90 1. of liquid nitrogen. Com- 
pressors 2RG3/350 and 1 VUV are used to compress the 
hydrogen. 


Table 1. Performance of the Liquefier 


Number 1 2 3 

Quantity of compressed 150 190 150 
hydrogen (m°(n.t.p.)/hr) 

Pressure of compressed 140 135 140 
hydrogen (atm) 

Pressure: 
In pumped nitrogen bath 140 140 140 
(mm Hg) 

In additional collector 0:3 0:5 2:5-3-0 
(atm) 

In the Dewar (atm) 0:3 0:5 0:3 

para-hydrogen in liquid 25 (normal) 25 (normal) 99 
produced (per cent) 

Output of liquefier (1./hr) 40-42 50-52 28-30 


We must point out that there is the possibility of boost- 
ing the liquefier by increasing the amount of compres- 
sed hydrogen to 260-280 m3(n.t.p.)/hr and by raising the 
pressure in the Dewar to 1-0 atm. The expected output is 
then ~ 70 1./hr of liquid hydrogen. 

The efficiency of the liquefier was also tested when 
filling the insulating space of the Dewar with carbon 
dioxide at 0-2 atm (after preliminary pumping out with a 
fore-vacuum pump). No change in the output was ob- 
served. 


The authors take this opportunity of expressing their 
thanks to V. A. Belushkin, L. B. Golovanov, L. P. 
Belonogova, and A. A. Belushkina, who took part in the 
construction of the liquefier; to R. A. Buyanov who made 
a study of the catalysis, and to K. A. Baicher, V. N. 
Dmitrievska, and other colleagues of the workshops and 
hydrogen station, by whose efforts the liquefier was built, 
set up, and tested. 
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IN an earlier publication! we described experiments from 


which we were able to deduce the electrical resistivity of 


the two phases of sodium at low temperatures. We were 
interested in the properties of the two phases of sodium 
partly because its low temperature phase provides the 
first known example of a hexagonal monovalent metal 
and partly because sodium has long been regarded froma 
theoretical point of view as a very good approximation to 
a free electron metal. Lithium also transforms at low 
temperatures by the martensitic mechanism, but in 
lithium there is an extra complication because the b.c.c. 
phase, which is stable at high temperatures, can transform 
to either of two faulted close-packed phases according to 
the conditions of the experiment.*~> Although lithium has 
only three electrons in the free atom and is thus funda- 
mentally a simpler element than sodium, its band struc- 
ture in the solid state is almost certainly more complicated 
than that of sodium and, at least in the body-centred phase, 
its Fermi surface is probably highly distorted.°* Lithium 
therefore provides an interesting contrast to sodium and 
it was our intention to determine the electrical resistivity 
of the three phases of lithium by experiments similar to 
those on sodium. For reasons which we explain below, 
we were not in fact able to do this in full detail, but we have 
established certain facts about the resistivity of the dif- 
ferent phases and these do indeed present a contrast to the 
corresponding data for sodium. 

It turns out that in general the electrical behaviour of 
lithium is more complex than that of sodium; for example, 
in its sensitivity to small amounts of impurity. We have 
therefore extended our measurements to include a range 
of dilute alloys of lithium with magnesium and by means 
of these measurements we have, we believe, been able to 
throw some light on the behaviour of pure lithium. We 
have reported elsewhere the results of our experiments on 
the electrical resistivity of lithium-6;? these experiments 
indicate that the effects of phase transformation are 
essentially similar tn the two isotopes. 
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The structure of the low temperature phases of lithium 


An important fact which emerged from the experiments 
of Barrett? and of Barrett and Trautz® is that the low 
temperature phase produced spontaneously by cooling a 
sample of lithium is predominantly close-packed 
hexagonal whereas that produced by deformation at low 
temperatures is predominantly close-packed cubic; in 
sodium only the close-packed hexagonai phase is formed. 
In writing of lithium, Barrett? says: ‘It would be appro- 
priate to describe the structure of the spontaneously 
formed low-temperature phase as a nearly random place- 
ment of each successive layer of atoms in the hollows of the 
layer just below, or alternatively as an f.c.c. structure with 
very imperfect stacking of (111) planes.’ Of the other close- 
packed phase he writes: ‘In the further transition to the 
f.c.c. structure the sequence [of close-packed layers] 
becomes A’B’C’A’'B'C’ . . . without change in the 
hexagonal atom pattern of each layer, and with relatively 
infrequent breaks in the sequence, a change brought 
about by deforming the metal at low temperatures.’ 


Details of the martensitic transformation in lithium 


We shall give here the basic facts about the martensitic 
transformation in pure lithium and in lithium—magnesium 
alloys; for a description of the general features of martensi- 
tic transformation we refer the reader to reference 1. 

Figure l(a) shows schematically the proportion of the 
low temperature phase as a function of temperature for a 
substance undergoing a martensitic transformation. 
Figure 1(b) shows the difference in phase composition as a 
function of temperature; it is the difference between 
curves (1) and (2) of Figure I(a) plotted against the 
temperature. This difference in phase composition shows 
itself as a temperature hysteresis in a number of physical 
properties (cf. Figure 3). The important temperatures 
from our point of view are M,, the temperature at which 
the transformation begins; M,, the temperature at which 
reversion begins; and 7g (or M,r- in some schemes of 
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notation) the temperature at which reversion is complete. 

Large grained specimens of pure lithium can have an 
M, as high as 80°K,’° and plastically deformed specimens 
can have M, as low as 67:°5°K.!! For pure lithium, how- 
ever, M, is usually between 70° and 75°K, and the speci- 
mens used in this work nearly all had values of about 72° K. 
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Figure I (a). Schematic diagram showing the fraction of the low 
temperature phase present at different temperatures in a substance 
undergoing a martensitic transformation 


Figure I (b). The difference between the phase composition of the 

substance at a given temperature on first cooling to that temperature 

and on re-warming to that temperature after cooling to very low 
temperatures 


The experiments on dilute alloys of magnesium in lithium 
show that M, increases by about 10° K for each atomic per 
cent of added magnesium.’ The reversion region in 
lithium (M, to 7g in Figure 1) is from 95° to 170°K, 
rapid reversion beginning at about 100°K and extending 
to about 120° K where the rate of reversion becomes small. 
In the alloys the bulk of the reversion appears to occur ata 
temperature which is higher than in pure lithium by about 
10°K per atomic per cent of added magnesium although 
the temperature of the end of the reversion (Tg) appears 
to be rather lower than in pure lithium, at about 160° K.!2 

The amount of the h.c.p. phase that can be produced by 
thermal cycling appears from X-ray studies on pure 
annealed samples of lithium to be about 50 per cent at 
about 60°K;*> the same appears to be true in lithium 
magnesium alloys. According to Hull and Rosenberg" 
about 90 per cent of a pure lithium specimen transforms 
spontaneously by 6°K; this would seem to be consistent 
with Barrett’s figure at 60°K if the transformation curve 
has the shape deduced by Martin'* from his specific heat 
measurements. This large amount transformed would 
indicate that the transformation is not readily stopped by 
strain. It is therefore not surprising that the amount trans- 
formed is not very sensitive to the physical state of the 
specimens used. Moreover, Martin’s experiments show 
that simple thermal cycling does not inhibit the trans- 
formation in lithium to any significant extent in contrast 
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to the behaviour of sodium.'!> Because we were able to 
inhibit the transformation in sodium by thermal cycling, 
we were able to study the resistivity of the two phases of 
sodium in considerable detail; the absence of this effect in 
lithium has prevented us from obtaining such detailed 
information about lithium. 

The effect of plastic deformation at low temperatures 
was studied by Barrett and Trautz? in pure lithium and in 
lithium—magnesium alloys. As already mentioned they 
found that plastic deformation produced an f.c.c. phase. 
Cold work at 78°K produced up to 50 per cent of this 
phase in pure lithium and higher proportions in dilute 
lithium—magnesium alloys. They found that the lower the 
temperature of deformation the greater the fraction of 
f.c.c. phase that was formed and although they made no 
deformation experiments below 78° K it seems reasonable 
from their work at higher temperatures to assume that at 
4°K almost the whole specimen can be converted by cold 
work to the f.c.c. phase (cf. Figure 8 of reference 5). 
In the present work we have deformed wires by stretching 
at 4°K and by squashing at 78°K. In the experiments at 
4°K we shall assume that in both pure lithium and the 
alloys we were able to convert the whole specimen to the 
f.c.c. phase; in the experiments at 78° K we shall see below 
that it is possible to estimate the amount of f.c.c. phase 
which has been produced by the deformation. All the 
estimates of the amounts transformed are approximate 
with an error of at least + 5 per cent, and more probably 
+10 per cent. 

We have considered the question of the distribution of 
the phases in our work on sodium.! In lithium we shall see 
that all the changes of resistivity on transforming are 
relatively small, and in this case we make no serious error 
in assuming that the various phases act as if they are 
electrically in series. 

Our description of the transformation in lithium has so 
far been based on fairly well established facts, but there 
are three other phenomena which may be important in 
these experiments but about which little can be said with 
certainty. 

Firstly, it is known from the experiments of Barrett and 
Clifton'® that in at least some lithium—magnesium alloys 
there can be thermoelastic changes of phase composition 
(i.e. a change of composition with temperature which 
does not show a temperature hysteresis). However, the 
alloy which they studied was relatively rich in magnesium 
and, for the same reasons that we have already given for 
sodium,! we shall assume that none of our specimens 
undergoes appreciable thermoelastic transformation. 
Secondly, we know that in some transformations of the 
martensitic kind the amount of transformation can be 
reduced by suitable deformation or heat treatment; this is 
known as stabilization. Although these effects are small in 
lithium, they must modify the transformation and rever- 
sion to some extent (cf. the reversion curve for lithium 
observed by Martin'*). Nevertheless, we shall ignore any 
effects of stabilization on our estimates of the amount of 
transformation (see, in particular, page 110). 
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Lastly, we have, in the outline above, taken the simple 
view that only the b.c.c. and h.c.p. phases are involved in 
the spontaneous (thermal) transformation and reversion. 
It appears that the thermodynamically stable phase is the 
f.c.c. one and it would therefore not be surprising to find 
that the f.c.c. phase can play some part even in the spon- 
taneous transformation. The X-ray results seem to show 
quite definitely that the f.c.c. phase is not produced during 
spontaneous transformation, but there is no firm X-ray 
evidence on whether or not it can exist as a transient phase 
during the reversion. The present experiments as well as 
those of Martin!* indicate that there may be some such 
effect in well annealed samples of pure lithium. 


The specimens 


In these experiments we studied samples from three 
batches of pure lithium which we shall refer to as A, B, and 
C; they were of increasing purity in that order, although 
B and C were rather similar. We also studied a series of 
dilute alloys of magnesium in lithium which were made 
from lithium of batch A with the appropriate amount of 
high purity magnesium added. A master alloy containing 
10 atomic per cent of magnesium was prepared by the 
method described by Barrett and Trautz,° and alloys of 
the required concentration were then made by diluting the 
master alloy with pure lithium. 

We shall distinguish between the different specimens of 
pure lithium by giving each one a number (which merely 
indicates the order in which the specimens were made) 
and a letter, A, B, or C, which indicates the batch from 
which the specimen was made. The alloys are distinguished 
by giving in each case the amount of magnesium in the 
alloy; thus the alloy containing 0-95 atomic per cent of 
magnesium we write as LiMg 0:95%. 

In order to compare the alloys with each other and with 
pure lithium it was necessary to determine their absolute 
resistivities at one temperature. This was done by extrud- 
ing a sample of the material into the form of a rod, about 
10 cm long and 0-5 cm in diameter. The dimensions and 
electrical resistance of the rod were then measured at 
room temperature. The mean diameter was measured with 
a micrometer, the average random error in these measure- 
ments being about 0:5 per cent and the estimated total 
error in the resistivity about 1 per cent. Values of the 
resistivity of pure lithium and of several alloys of known 
composition are given in Table 1 together with the cor- 
responding residual resistivity, p9, measured at 4:2°K. 
Details of all the pure lithium samples are given in Table 2. 


Experimental results—the h.c.p. phase 


To study the electrical resistance of the spontaneously 
formed h.c.p. phase, we used two methods which we shall 
call ‘the variable temperature method’ and ‘the reference 
temperature method’. In the first, we cool the specimen to 
a low temperature (4° or 10°K), measuring the electrical 
resistivity at about 10° intervals during the cooling process 
and then again at the same temperatures during the warm- 
ing process. In this way we find out how the resistance of 
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the specimen at various temperatures varies with the 
proportion of the two phases. Unfortunately this propor- 
tion is also a function of temperature and the two effects 
are hard to separate. In the second method, the resistance 
of the specimen is always measured at a fixed reference 
temperature above the M, temperature and below that 
at which reversion begins. Between measurements the 
specimen is cooled to successively lower temperatures to 
study the effect on the resistance at the reference tempera- 
ture of different proportions of the two phases. Because 
the reference temperature is below the reversion tempera- 
ture the phase composition is not appreciably altered in 
warming the specimen back to the reference temperature 
and so the phase composition of the specimen corresponds 
to that of the lowest temperature to which it has been 


Table 1. The Resistivity of Lithium—Magnesium Alloys 


Poo5 Pot Magnesium in lithiumt 
(Qem) (Qem) (atomic per cent) 
9-52 x 10-8 13:Srt On” (0), pure lithium, batch B 
9-58 59-5 0:05+0-01 
— 78 0:07+0-01 
9-78 250 0:21+0-01 
10:27 672 0:-47+0-015 
10-68 1036 0-84+0-015 
11-07 1300 0-95+0-015 


+ These results refer to two-phase mixtures (predominantly h.c.p.). 
} We are indebted to Dr. Inman of the Department of Mines, Ottawa, for the 
flame-photometric analysis of these alloys. 


Table 2. Summary of the Results of the Experiments on Pure 


Lithium 
| Percentage 
Di Bed: Residual change at 
Speciment AINE CEL HS COUT resistivity | 80°K due 
(mm) of Q 
specimens (Stem) abe 2 
formation 
Li (7A) 0-5 Nonet | 25-6~x 10-° —0-85 
Li (16A) 0:5 None 30 —0-28 
Li (8B) 0:5 None 13-8 +0-:76 
Li (12C) 0:5 None 10-6 + 0:66 
Li (13C) 0:5 None 10-6 +0-70 
Li (15C) 0:5 None 10-9 + 0-68 
Li (19C) 0:5 None 10-6 + 0-68 
Li (18C) 3 20 hr 8-9 —0-27 
at 150°C 
Same specimen 6 hr Shp 0-00 
after twice at 50°C 
cooling to 4°K 
Li (17C) 5 24 hr 8-2 0+0-2§ 
at 150°C 


+ The letter after the specimen reference number indicates the batch from which 
the specimen was made. Batches A and C were obtained from the Lithium 
Corporation of America. Batch B is of unknown origin. 

{ This means that, after extrusion, the wires were kept at room temperature for an 
hour or two before measurements began. 

§ The specimen had such a large diameter that very accurate measurement of the 
resistance was not possible. 
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cooled. Then the specimen is warmed to successively 
higher temperatures between readings to study the rever- 
sion process. 


0) 50 100 150 200 
Temperature ——> (°K) 


O LiMg0:21% @ LiMg0:95% 
Figure 2. The change in resistivity of two lithium—magnesium alloys 


at 80:6° K as a function of the temperature to which the specimens were 
successively cooled (or heated) between readings 


(1) The alloys 


We consider first the results for the alloys since these 
are simpler. 

The results of the reference temperature method applied 
to LiMg 0-21°% and LiMg 0-95°% are shown in Figure 2. 
If, as we might expect, the changes in the resistance at the 
reference temperature (~ 80° K in these experiments) are 
proportional to the fraction of the h.c.p. phase present, 
then these curves should have the same shape as the trans- 
formation curves of the two specimens (cf. Figure 1(a)). 
The shape of these curves does indeed agree fairly 
well with the transformation curves deduced by other 
methods;*!* we may therefore feel confident that we 
understand what is happening at this temperature. Two 
other features of the curves are noteworthy: first, that the 
resistivity of the transformed material is lower than that 
of the untransformed and, second, that the magnitude of 
the change is about the same in the two alloys. Expressed 
as a fraction of the ideal resistivity at this temperature, the 
three alloys (LiMg 0-21%, 0-47%, and 0-95%) all undergo 
a decrease of resistivity of between 5 and 8 per cent due to 
the transformation. 

We consider now the results of the variable temperature 
method. Figure 3 shows the percentage difference between 
the resistance of the alloy specimens measured with 
decreasing temperature and that measured with increasing 
temperature, this difference being plotted as a function of 
the temperature. If the resistance due to impurities and 
that due to phonons both changed in the same way during 
the transformation and if Matthiessen’s rule were applic- 
able, the results of the variable temperature method would 
be closely related to those of the reference temperature 
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method. The curve resulting from the former method 
would then roughly resemble that of Figure 1(b). A com- 
parison of Figure 2 with the corresponding two curves in 
Figure 3 shows that above about 70°K the two different 
methods both show that the h.c.p. phase has a lower 
resistance than that of the b.c.c. phase. At lower tempera- 
tures, however, where the residual resistivity dominates, 
Figure 3 shows that the transformation now increases the 
resistance of the specimen, in complete contrast to the 
behaviour at 80°K shown in Figure 2. As we shall show 
below, this difference is attributable to departures from 
Matthiessen’s rule. 

It is possible to get a rough estimate of the magnitude of 
the change in residual resistance produced by the spon- 
taneous transformation. We can use the information in 
Figure 2 as a guide to the amount of transformation at 
different temperatures. (If we use the transformation 
curves derived from specific heat measurements!* we get 
similar results.) We see from Figure 2 that at 40°K the 
difference in composition on cooling and warming is 
rather less than 10 per cent of the total change which 
occurs. If the total amount spontaneously transformed at 


Sats 


Percentage change in resistance 


0) 50 100 150 200 


Temperature ——> (°K) 


A LiMg 0-21 WA O LiMg0-47% V LiMg 095% 
© Pure lithium specimen 7A [J Pure lithium specimen 8B 
Figure 3. The percentage difference between the resistance of specimens 
measured with decreasing and increasing temperatures plotted as a 
function of temperature (the negative sign indicates a lower resistance 
after transformation) 


CRYOGENICS : DECEMBER 1961 


4°K is about 90 per cent, the change in composition at 
40° K (before and after cooling to 4:2°K) is thus about 10 
per cent. According to Figure 3, this difference in composi- 
tion increases the resistance of LiMg 0-95°%% by about 
0-5 per cent, from which we deduce that a complete trans- 
formation would increase the resistance by 5 per cent. 
Since at 40°K most of the resistance is residual we deduce 
that the residual resistance of the spontaneously formed 
h.c.p. phase is at least 5 per cent higher than that of the 
b.c.c. phase. A similar result follows for the LiMg 0-21% 
alloy. Both these results tend to underestimate the change 
since the effect of the phonon scattering, which is in the 
opposite direction to that of the residual scattering, has 
not been taken into account. Our conclusion therefore is 
that the residual resistivity in the h.c.p. phase of the 
alloys is 5-10 per cent higher than that in the b.c.c. phase. 


(2) Pure Lithium 


Our results for several different specimens of pure 
lithium have been collected in Table 2. The last column 
gives the percentage change in resistance measured at 80°K 
caused by the spontaneous transformation (to transform 
the specimens, they were cooled to 4-2°K and brought 
back to 80°K). It is seen that all the changes are quite 
small. Certain other general features are clear: (1) the 
specimens from batch A, which are less pure than the 
others, show a decrease in resistance at 80°K on trans- 
forming; (2) specimens of the same size and heat treat- 
ment, but from the purer batches B and C, consistently 
show an increase of resistance-on transforming; (3) the 
purest material (batch C) can be made to show a decrease 
of resistance on transforming if the specimens are big 
enough and are well annealed (e.g. specimen 18C; see 
also Table 3). Results of measurements by the variable 
temperature method on two pure lithium specimens have 
also been included in Figure 3. 5 


Table 3. Resistance (in Arbitrary Units) of Lithium Specimen 18C 
after Annealing at ~50°C for 6 hr 


Temperatures at which resistance 
was measured (° K) 
34:2 S501 63:8 70:3 80:6 


Temperatures to 
which the specimen 
was cooled (° K) 42 


80-6 a — — — — 16035 
70°3 = == — — 11015 15977 
63°8 — —_ — 8266 10996 15977 
50:1 — — 3824 8281 11013 15983 
34:2 — 1041 3833 8297 11038 16006 
42 143, — — — 11043¢ 16037f 


+ Subsequent run after re-anneal. 


Before going into further details we will give what we 
believe are the reasons for this behaviour. At temperatures 
at which the resistivity is predominantly due to phonons, 
the primary effect of the transformation, as in the alloys, 
is to decrease the resistivity of the specimens. However, for 
reasons which we shall discuss later, the magnitude of this 
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decrease is much smaller in the pure specimens than in the 
alloys. In the purest specimens the decrease may be so 
small that if the transformation itself generates a large 
number of new scattering centres, thereby substantially 
increasing the residual resistivity of the specimen, this 
secondary effect can outweigh the primary effect and a net 
increase of resistivity can result. We expect that, in a large 
well-annealed specimen, the regions which transform (at 
least in the early stages of the transformation) will be com- 
paratively big and free from defects in contrast to those ina 
small-grained specimen of small diameter. We thus expect 
the transformation to increase the residual resistivity of 
the small unannealed specimens much more than those 
which are large and well annealed. This seems to be borne 
out by the results in Table 2. 

Another consequence of this explanation is as follows. 
From the nature of a martensitic transformation it is 
known that in general the specimens tend to become more 
strained as the transformation proceeds and that the 
regions which first transform are usually the largest and 
the most strain-free. This would mean that the specimen 
resistance as measured at a given temperature would tend 
to fall during the early stages of the transformation and 
tend to increase during the later stages. This kind of 
behaviour was indeed observed in an experiment on 
specimen 18C. In this experiment the specimen was cooled 
to successively lower temperatures below 80° K and after 
each cooling the resistance was measured at a series of 
reference temperatures up to and including 80°K. The 
temperatures and the values of resistance are given in 
Table 3. It is seen, for example, that as the proportion of 
the low temperature phase increases, the resistance at 
80-6°K first falls and then rises. The resistance of the 
specimen at 70-3° K does likewise, but at 63-8°K, at which 
temperature only the later parts of the transformation 
take place, no decrease of resistance is seen. 

Details of the reversion of specimen 18C are given in 
Figure 4 while in Figure 5 are shown the complete con- 
version-reversion cycles of specimens 7A and 8B, all 
obtained by the reference temperature method. (These 
figures should be compared with Figure 2 for the alloys.) 
During the reversion cycle all the pure specimens show a 
‘hump’ just above 110°K. This is almost certainly due to 
an increase in the residual resistivity caused by new 
strains and defects generated in the specimen by the 
reversion processes. (This hump can also just be seen in the 
results of the variable temperature method although the 
effect is then found to be much smaller.) Ultrasonic 
experiments on samples of pure lithium!! show that at 
this temperature during the reversion process there is a 
marked increase in the attentuation of the ultrasonic 
waves. The large peak which occurs at 160° K in the rever- 
sion curve of specimen 18C and to a lesser extent in that of 
specimen 7B is difficult to account for; since it occurs at a 
temperature at which almost all the low temperature phase 
is known to have disappeared it must imply that the 
resistivity of the body-centred phase has been substantially 
increased in the transformation and reversion process. 
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(3) The change of shape of specimens on transformation 


The experiments of Basinski and Verdini'’ have shown 
that the change of specific volume on phase transformation 
is small. On the other hand, we know that the shear 
strain associated with the transformation is fairly large (cf. 
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Figure 4. The change in resistivity of pure lithium specimen 18C 
measured at 80° K as a function of the temperature to which the 
specimen was successively heated between readings 


the discussion in the paper by Basinski and Verdini). We 
have therefore made some measurements to see how the 
change of shape of the wires we have used can affect the 
measured values of resistance. We measured directly the 
change in length of some thin wires which were prepared 
in the same way as the wires used in the resistance work and 
it appears that the change of shape in transforming can 
be sufficient to increase the resistance by about 0-5 per 
cent quite apart from any changes of resistivity which 
may occur. The shape factor, moreover, does not change 
in proportion to the phase composition of the wire, and in 
some cases about half of the change of shape factor per- 
sists after the specimen has completely reverted. 

In reducing the resistance results to values of resistivity 
we have assumed that the shape factor of our wires does 
not change with transformation. Thus this effect of the 
change of shape factor will be to give an apparent increase 
of resistivity. This effect shows itself to a certain extent in 
Figures 3, 4, and 7, for instance, where there are indications 
of a change of shape factor of about 0-1 per cent as a 
result of the transformation. This is smaller than the 
experiment on the change of length would lead us to 
expect, and it indicates that the change of shape factors are 
only important in our thin wires at the 0-2 per cent level. 
In the thicker rods, where the preferred orientation of the 
crystal grains is less, the effect should be less serious. This is 
probably a partial explanation of the change in sign of the 
results in the last column of Table 2 as the specimen size 
is changed. 
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This uncertainty about the shape factor means that the 
changes of resistivity on transformation are uncertain to 
about 0:2 per cent. This does not invalidate the conclusions 
we shall draw here, but it does mean that in a detailed inter- 
pretation of the reversion curves (cf. Figure 5) this extra 
effect must also be considered. 


(4) Departures from Matthiessen’s rule 


From the measurements on the dilute lithium—magnes- 
ium alloys it soon became evident that Matthiessen’s rule 
was not obeyed. In order to show this, we have subtracted 
from the measured resistivities of the alloys at each 
temperature the corresponding ideal resistivity as deduced 
from measurements on samples of the purest lithium. The 
ideal resistivity of pure lithium may itself be altered by the 
transformation but this uncertainty is small enough to be 
neglected for our present purpose. 

In this way we have found out how the apparent residual 
resistivity, po, varies with temperature for the three 
alloys LiMg 0:21%, 0:-47%, and 0-95%. The results are 
shown in Figure 6 as a function of both decreasing and 
increasing temperature. If Matthiessen’s rule were 
strictly obeyed and if there were no transformation, po 
would be a constant independent of temperature. 

It is seen that below 100°K, po falls markedly with 
falling temperature. This cannot be ascribed to the change 
of residual resistivity on transformation, since as we have 
seen from the variable temperature method, the trans- 
formation causes an increase in the resistivity at low 
temperatures. In other words, if the transformation did 
not occur the low temperature resistivity of the alloys 
would be even smaller thereby further decreasing the 
value of po at low temperatures. It is therefore clear that 


0 0) 100 ISO 200 250 300 
Temperature =——> (°K) 
—@— Pure lithium specimen 7A measured at 74:8° K 


--O©O-— Pure lithium specimen 7A measured at 90-4°K 
—-[-— Pure lithium specimen 8B measured at 75-5° K 


Figure 5. The change in resistivity of pure lithium specimens at a fixed 
reference temperature as a function of the temperature to which the 
Specimens were successively cooled (or heated) between readings 


CRYOGENICS - DECEMBER 1961 


in the b.c.c. phase there are large departures from Mat- 
thiessen’s rule. 

If at 4°K the specimens are almost completely trans- 
formed, the py curve measured during the warming pro- 
cess is, up to about 90°K, essentially that of the h.c.p. 
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Figure 6. The apparent residual resistivity of lithium—magnesium 

alloys as a function of temperature, illustrating departures from 

Matthiessen’s rule. The open symbols refer to measurements made 

with decreasing temperature and the full symbols to measurements 
made with increasing temperature 


phase. It is thus evident that departures from Matthies- 
sen’s rule must occur in this phase too, although they are 
smaller than those found in the body-centred phase. 

We have also compared the resistivities of specimens 
from two batches (A and C) of pure lithium and we find 
that marked departures from Matthiessen’s rule occur 
here also as illustrated in Figure 7. In this we plot Ap, the 
difference in the resistivities of the two specimens as a 
function of decreasing and increasing temperature. We 
have indicated schematically how Ap for the body- 
centred phase might be expected to vary with temperature 
if it could be measured down to the lowest temperature 
(dashed line in Figure 7). 


Experimental results—h.c.p. phase; summary 


The main features which emerge from these experiments 
may be summarized as follows: 


(1) The spontaneous transformation from the b.c.c. to 
the h.c.p. phase increases the resistivity of the alloys at 
low temperatures. 

(2) The spontaneous transformation always decreases 
the resistance of the alloys at 80°K. 

(3) In pure samples of lithium the transformation can 
either increase or decrease the resistance at 80°K but all 
the changes are quite small. In a thick well annealed 
sample of the purest lithium the decrease in resistance at 
80° K is about 0-25 per cent. 

(4) There are large departures from Matthiessen’s rule 
in both phases of the alloys and at least in the b.c.c. phase 
of quite pure lithium. 
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Experimental results—f.c.c. phase 


(1) The straining experiments at 4:2°K 


In order to study the residual resistivity of the f.c.c. 
phase, bare wires of pure lithium and of the lithium— 
magnesium alloys were deformed by stretching them in 
liquid helium. The resistance was measured as a function 
of strain and the results are shown in Figures 8(a) and 
8(b). These show the resistivity of the specimens as a 
function of In(///o), where / is the instantaneous and Jy 
the initial length of the specimen. In these curves we see 
that there are two regions of distinctly different character: 
first, the initial non-linear change of resistivity with strain 
and, second, the almost linear region which occurs at 
somewhat higher strains. 

By comparison with the results of straining experiments 
on sodium! and on potassium!® we attribute the almost 
linear region to the increase with strain of the number of 
dislocations in the specimens, whereas the curved region 
at low strains we attribute to the martensitic transforma- 
tion, probably a transformation from the h.c.p. to the 
f.c.c. phase. This curved region is thus the part which 
interests us here and we have estimated the change in 
residual resistivity associated with this transformation, 
Ap’, by extrapolating the linear part of the curve to zero 
strain as indicated in the figures. These values are listed in 
Table 4. 

From the data shown in Figure 8 and in Table 4, we see 
that the transformation increases the residual resistivity ; 
the values of Ap¢'/po for the alloys are all about 0-06 while 
the corresponding value for pure lithium is considerably 
greater. We discuss possible interpretations of these 
results below. 


0 50 100 oe os 700 
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Figure 7. The difference in resistivity of two pure samples of lithium 

as a function of temperature. The sharp cusp at about 70°K in the 

falling temperature curve shows the onset of the transformation. The 

dashed line indicates schematically how the difference in resistivity of 

the pure body-centred phases might change with temperature below the 
transformation temperature 
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Table 4. Summary of the Straining Experiments 


assumptions we make are as follows: (i) that the de- 
formation produces a fraction f of the f.c.c. phase only 
(we checked by measurements on the alloys at 90° K that 
no spontaneous transformation occurred before deforma- 
tion); (ii) that during the subsequent cooling to 4°K, the 
remaining fraction (1 —f) of the b.c.c. phase is converted 
spontaneously to the h.c.p. phase while the f.c.c. phase 
remains unchanged; (iii) that in the final stage of the 
experiment, in which we start with the same specimens 


Approximate Ao” x 
atomic percentage of p3(Qem) mal 703 (Qem) 
magnesium in lithium Po . 
(0) 15-4x10-* 0-11 235+5x 10° 
(0) 14-6 0-17 _— 
0-07 82:8 0:06 335+10 
0:35 466 0:05 505 +50 
0-73 969 0-065 515+50 
0-95 1290 0-058 550+ 100 


(2) Deformation experiments at about 80°K 


In order to make a complete study of the temperature 
dependence of the resistivity of the f.c.c. phase it would be 
necessary to use a deformation cryostat capable of operat- 
ing over a wide temperature range. It is, however, possible 
to obtain a certain amount of information about the f.c.c. 
phase from experiments using suitable fixed temperature 
baths and this is what we have done. The experiments 
consisted of first squashing a piece of wire (lithium or 
alloy) under liquid oxygen or nitrogen using a pair of pre- 
cooled pliers; in this process the amount of deformation 
corresponded to an average of about 50 per cent reduction 
in thickness. The wire was then immersed in a series of 
different temperature baths (liquid helium, nitrogen, oxy- 
gen, acetone, and solid carbon dioxide) and the effect on 
its resistance of various thermal treatments was measured. 
We also made a number of experiments to test the reliabil- 
ity of these results and to make sure that the transient 
heating effects when the specimen was moved from one 
bath to another did not affect the results. 

The sequence of the measurements and the assumed 
proportions of the different phases at each stage are 
shown in Table 5. As the table shows, the important 


Table 5. Sequence of Measurements in the 
Deformation Experiments 


Composition 
b.c.c he.p. | fiec 
Squash at 77°K (or just above M,) ! 
(1) Measure at 77°K, 90°K (1-f) 0 i 
(2) Cool and measure at4:2°K . 0O (1-f) if 
(3) Remeasure at 77°K,90°K =§ O (1-S) if 
Make specimen revert at 195°K (i.e. | 
do not anneal) | 
(4) Measure at 77°K, 90°K 1 0 0 
(5) Cool and measure at 4:2°K 0 1 0 
(6) Remeasure at 77°K,90°K 0 1 0 
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Figure 8(a). The resistivity at 4:2°K versus strain for specimens of 
pure lithium (B) (scale A) and for LiMg 0-07% (scale B) 


Scale B Scale C 
tl. eee te eee < 
o LiMg 0:95%| L- B as 
a | Scale C4 ae = ~ 

Veen ee +1100, 7 4400 
ae LiMg 0:73% x10 x10 
| Fat ee Scale B | 
550} (== \ | 
x10 Bae 7 
f foi inh 
Gin 4 i | vs a < < 
iy, ; | 1,000 1,300 
y Jf Lima 0:35% 
a Scale A 


0.05 


0:10 0: 


Strain Un (2/1) 


1,200 


\S 


ae 


Figure 8(b). The resistivity of three lithium-ma i : 
, enesium alloys: 
LiMg 0:35 % (scale A); LiMg 0-73 % (scale B); LiMg 0:95 % (scale C) 
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untransformed at 77° or 90°K, cooling to 4°K converts 
them virtually completely to the h.c.p. phase. 

Consider first the total resistance of the specimen at 
either 77° or 90° K. The effect of complete transformation 
from b.c.c. to h.c.p. is given by the difference between 
measurements (4) and (6) of Table 5. On the other hand, 
the effect of transforming a fraction (1 —/f) of the specimen 
from b.c.c. to h.c.p. is given by the difference between (1) 
and (3). This therefore enables us to estimate f, the amount 
of f.c.c. phase formed by the deformation. If we now 
assume that strain contributions to the resistivity are the 
same at (1) and at (4), we can deduce from the difference 
between (1) and (4) the value of the resistance of a speci- 
men containing only the f.c.c. phase. 

This assumption about the strain contribution implies 
that by warming the specimens to 195°K we have caused 
them to revert to the b.c.c. phase without appreciably 
annealing them. This can be tested from the measurements 
at 4°K. Since we know the value of f at 4°K from the 
measurements at 77° or 90°K, we can use the result of 
section (1) above (page 109) to calculate what difference 
we should expect to find between the measurements (2) 
and (5), assuming that the same scattering centres exist in 
both states of the specimen. According to our findings in 
section (1), the intrinsic residual resistivity of the f.c.c. 
phase is about 6 per cent greater than that of the h.c.p. 
phase. Using this information we find that in our present 
experiments the observed difference in resistance between 
states (2) and (5) is close to that calculated on the assump- 
tion that the strain contribution to the resistivity is about 
the same for the two measurements. Our assumption is 
thus justified. 

The results for f and for the resistance of each of the 
separate phases are given in Table 6 for two lithium— 
magnesium alloys; the values of fso found appear quite 
reasonable? and the two values for each alloy show that the 
measurements at the nitrogen and oxygen points are 
consistent. The same analysis applied to the measure- 
ments on pure lithium gave results with internal incon- 
sistencies but in view of the results of section (1) it does not 
seem unreasonable to suppose that these may be due to 
the less predictable effects of strain in pure lithium and 
are not due to any fundamental error in the method. 


Discussion 


The following is an outline of the argument which we 
present in detail in the subsequent sections. In section (1) 
we conclude from the results of various experiments on 
lithium and dilute lithium—magnesium alloys (Knight 
shift, thermoelectric power, electrical resistivity, specific 
heats) that the properties of the Fermi surface at least in 
the b.c.c. and h.c.p. phases are not much changed by 
alloying with up to 1 per cent of magnesium. In section 
(2) we deduce from the results of the straining experi- 
ments at 4°K that there is an intrinsic difference between 
the electron properties of the h.c.p. and f.c.c. phases. In 
section (3) we argue that there is also a difference between 
the electron properties of the b.c.c. and h.c.p. phases and 
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that this may explain why the b.c.c. and h.c.p. phases of 
pure lithium have almost the same resistivity at 80° K, in 
spite of their different Debye temperatures. In section (4) 
we observe that, because departures from Matthiessen’s 
rule are smaller in the h.c.p. phase than in the b.c.c. phase, 
this difference between the phases probably accounts for 
the much larger changes in resistivity that are observed at 
80° K when the alloys transform than are observed when 
pure lithium transforms. In section (5) we consider what 
we can deduce about the resistivity of the f.c.c. phases of 
the lithium—magnesium alloys at 80°K. 


(1) Change of the Fermi surface on alloying 


Hughes!’ has made measurements of the Knight shift in 
lithium and lithium—magnesium alloys from which it 
appears that the density of states in lithium is not greatly 
changed by alloying with magnesium. These measure- 
ments, made on the body-centred phase, show that the 
Knight shift decreases initially by about 1 per cent per 
atomic per cent of magnesium. The change in density of 
states, however, is thought to be considerably less than 
this.!? We shall therefore assume that in the dilute alloys 
with which we are concerned the addition of magnesium 
does not significantly alter the Fermi surface in the body- 
centred phase. 

Martin”? has measured the electronic specific heat of 
both lithium and some lithium—magnesium alloys. His 
measurements were made on spontaneously transformed 
specimens at about 1°K which would thus be predomin- 
antly in the h.c.p. phase. His results show that the pure 
lithium specimen and the LiMg 1% alloy had the same 
electronic specific heat to within the experimental error 
of about 1 per cent. Moreover, experiments on the 
thermoelectric power below 1°K of lithium and its dilute 
alloys with magnesium”! show no sign that the Fermi 
surface is changed by alloying at least up to 1 per cent of 
magnesium. In the discussion which follows we shall 
therefore assume that in none of the phases of lithium is 
the Fermi surface appreciably altered by the alloying with 
magnesium. 


(2) The change of residual resistivity in the strain-induced 
transformation 


The change of residual resistivity on transforming, 
Apo’, aS measured by the straining experiments, could 
arise from an intrinsic difference in the resistivity of the 
two phases (i.e. from a difference in electron properties 
in the two phases) or from changes in the number of im- 
perfections capable of scattering the electrons or from a 
combination of both of these. The number of scattering 
centres does indeed change during the transformation; 
according to Barrett the large number of stacking faults 
in the h.c.p. phase is greatly diminished when the f.c.c. 
phase is formed (see page 103). On this basis we might 
expect the resistivity to diminish during the transforma- 
tion from h.c.p. to f.c.c. However, as we have noted else- 
where, the scattering from stacking faults in pure lithium 
is almost certainly very small;?? this would probably be 
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true also of the alloys so that we shall ignore the effect of 
stacking faults here. 

The values of Ap3"/po (see Table 4) for the alloys are all 
fairly constant and equal to about 0-06, even though the 
magnitude of the residual resistivity itself changes: by a 
factor of 15. Now we believe that the effect of tensile 
strain in transforming the specimens at 4°K is about the 
same in all the alloys, viz.: that we start with specimens 
predominantly in the h.c.p. phase and that strain converts 
them to predominantly the f.c.c. phase (cf. page 103). 
Thus the constancy of Ap@'/py would seem to imply that 
the change is due to the different properties of the electrons 
in the two phases. 

If we now consider the behaviour of pure lithium, we 
find that the value of Ap@'/po is considerably higher than 
for the alloys. On the basis of our assumption discussed 
above, that the Fermi surface is not appreciably changed 
by alloying, we interpret this as follows. We assume that 
in pure lithium, as in the alloys, there is a contribution to 
Ap9'/po of 0-06 from the change in electron properties, and 
that the remainder arises mainly from scattering by 
additional physical defects. (This scattering is presumably 
also present in the alloys but, assuming that it has approxi- 
mately the same magnitude as in pure lithium, its contri- 
bution would be negligible in the more concentrated 
alloys. Even in the most dilute alloy it would contribute 
only about 25 per cent to Apg' and this could not be 
recognized with any certainty.) If this change of 0-06 
corresponds to a complete change from h.c.p. to f.c.c. 
then the ratio of the residual resistivities in the two phases 
polf.c.c.)/po(h.c.p.) is equal to about 1-06. 


(3) The change in resistivity in the spontaneous transforma- 
tion 

The residual resistivity. The spontaneous transformation 
increases the resistivity of the alloys at low temperatures. 
Since this occurs at temperatures where the scattering is 
almost all residual scattering, we conclude that the residual 
resistivity is higher in the h.c.p. phase than in the b.c.c. 
phase. This could arise if (1) there were an increase in the 
number of scattering centres in going from b.c.c. to 
h.c.p., and/or (2) if there were a change in the electron 
properties in going from one phase to the other. It is not 
possible to be certain which of these predominates but 
it does seem unlikely that additional physical defects 
could produce an increase in residual resistivity as large 
as 7x 10~® cm (estimated from the data on the most 
concentrated alloy), which is five times larger than the 
total residual resistivity of our purest specimens. If it can 
be argued that the scattering from stacking faults, of 
which there is an abundance in the h.c.p. phase, can some- 
how be enhanced by the presence of impurity ions then 
this increase might be attributed to the formation of 
stacking faults. It has to be remembered, however, that 
the stacking faults are formed at a low temperature where 
appreciable diffusion is not possible (hence the martensitic 
nature of the transformation). We shall therefore tenta- 
tively assume that a change in electron properties is 
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responsible for the observed increase in residual resistivity 
in going from the b.c.c. to the h.c.p. phase of the alloy. 
(Such a change in electron properties could not be a 
simple consequence of the volume change due to the trans- 
formation since this is known to be very small.'”) 

We have, in addition, independent direct evidence for a 
change in the band structure of pure lithium associated 
with the spontaneous transformation. Recent measure- 
ments of the magnetic susceptibility of lithium in the 
neighbourhood of 80°K both before and after the spon- 
taneous transformation show that the susceptibility 
(which is that of a paramagnetic) is greater in the h.c.p. 
phase than in the b.c.c. phase.?* 

The change in resistivity of pure lithium on transforming 
at 80° K. It is clear from our experiments that in pure 
lithium the changes in resistivity at 80° K due to the spon- 
taneous transformation are very small, probably less 
than 0:5 per cent. We know, however, from Martin’s 
calorimetric measurements‘ that the Debye temperature 
6, is about 2:5 per cent higher in the hexagonal close- 
packed phase of lithium than in the body-centred phase 
and we should therefore expect to see a correspond- 
ing change in the ideal resistivity, as indeed we do in 
sodium.!!5 

We should expect, judging from the behaviour of 
sodium where our experiments show that the electron 
properties are not changed by the transformation, that at 
80°K the resistivity of the hexagonal close-packed phase 
of lithium would be roughly 8 per cent lower than that of 
the body-centred phase because of the difference in 
6 values. The most probable reason why no such change is 
observed is that the change in electron properties, which 
in the preceding. sub-section we deduced from the low 
temperature resistivity changes in the alloys and which 
would tend to increase the resistivity by 5-10 per cent, 
almost entirely offsets the effect at 80°K of the change in 
the value of 0. 


(4) The change in resistivity of the alloys on spontaneous 
transformation at 80° K 


The change in 6p due to the transformation is about the 
same in the alloys as in pure lithium!? and the change in 
electron properties is also assumed to be the same (see 
sub-section (2) above). Wethereforeexpect that the change 
in ‘ideal’ resistivity (i.e. that part of the resistivity which 
arises from the phonon scattering) due to the transforma- 
tion will be about as small in the alloys as it is in pure 
lithium. If, therefore, Matthiessen’s rule were obeyed in 
both phases the change at 80°K would just reflect any 
change in the residual scattering. In fact the changes 
observed at 80°K are quite different and even in the 
opposite sense to those at 4° K (see page 106). We attribute 
this to departures from Matthiessen’s rule which are 
greater in the b.c.c. phase than in the close-packed phase 
(cf. Figure 6). Departures from Matthiessen’s rule not 
unlike those observed here have been found in dilute gold— 
silver and silver-gold alloys by Krautz and Schultz. 

Our observation that the change of absolute resistivity 
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on transforming was about the same in all the alloys 
appeared to indicate that the changes were chiefly in the 
ideal resistivity. We now believe that this is a fortuitous 
appearance which arises in the alloys we have studied 
only because the magnitude of the departure from Mat- 
thiessen’s rule is about the same for all these alloys. 

In Figure 9, we show schematically what we believe to 
be the behaviour of the apparent residual resistivity of the 
two phases as a function of temperature. It is seen (Figure 
9(a)) that at the lowest temperature po(b.c.c.) is less than 
po(h.c.p.). At high temperatures po(b.c.c.) is considerably 
larger than po(h.c.p.) because of the larger departures 
from Matthiessen’s rule. In Figure 9(b) we indicate how a 
specimen following the transformation and reversion 
curves of Figure 9(c) would behave. 


(5) The change of resistivity of the alloys on strain- 
induced transformation at 80° K 


The results given in Table 6 show that the resistivity of 
the f.c.c. phase in the alloys is slightly higher than that of 
the h.c.p. phase. It is impossible to interpret these without 
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Figure 9. Schematic illustration of the variation of the apparent 

residual resistivity of the alloys with temperature in the two phases. 

(a) The apparent residual resistivity of the two pure phases. (b) The 

apparent residual resistivity (heavy line) which a specimen might show 

in following the transformation-reversion curve in (c) (cf. Figure 6). 

(c) The fraction of the close-packed phase present as a function of 
temperature on cooling and warming 
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knowing the resistivity of pure lithium in the f.c.c. phase 
and this we have not been able to find from our experi- 
ments. However, we have never observed a large change 
of resistivity on deforming pure lithium at 80°K and it 


Table 6. Summary of the Deformation Experiments at 
About 80°K 


LiMg 0-21% LiMg 0-95% 


Nz point O, point | Nz point Oz point 


Resistance (milliohms) 


b.c.c. phase 10-33 ist 22:04 25-56 
h.c.p. phase 9:93 12-97 21:66 Dail 
f.c.c. phase 10:06 13-14 21-78 25:19 
f, amount of f.c.c. 50:30 0:30 | 0-74 0-76 
phase produced by 
deformation 


would not be inconsistent to suppose that the ideal 
resistivity of the f.c.c. phase is quite similar to that of the 
b.c.c. and h.c.p. phases at 80°K. If this is the case then as 
with the other two phases, the f.c.c. phase of the alloys 
also shows large deviations from Matthiessen’s rule. 


Summary 


The changes in electrical resistivity of sodium and 


‘lithium which arise from the low temperature trans- 


formation provide an interesting contrast. At 44°K, the 
resistivity of the body-centred phase of sodium exceeds 
that of the hexagonal close-packed phase by about 15 per 
cent. At 80° K, which is the corresponding temperature for 
lithium, the difference between the ideal resistivity of the 
same two phases of lithium is probably less than 0:5 per 
cent. On the other hand the transformation has no 
intrinsic effect on the residual resistivity of sodium (within 
the experimental error of + 1 percent), whereas in lithium 
we estimate that the resistivity of the hexagonal close- 
packed phase exceeds that of the body-centred phase by 
about 10 per cent. Finally, there is in lithium the f.c.c. 
phase, which has no counterpart in sodium. Its residual 
resistivity is about 10 per cent higher again than that of 
the h.c.p. phase. 
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A Liquid Helium Level Indicator 


(AN EK ° Rovinskit All-Union Electrotechnical Institute, U.S.S.R. 


Received 18 July 1961+ 


THE idea of using the phenomenon of the transition of pure 
tantalum into the superconducting state at 4:3°K for 
measuring the level of liquid helium was proposed by 
Feldmeier and Serin and first applied by Rasor? in an 
apparatus for continuous measurement of the quantity of 
liquid in a closed vessel. In a liquid level indicator with a 
feeler in the form of a probe it is desirable to reduce the 
dimensions and mass of the immersed tantalum element, 
so that it is advantageous to change the electrical circuit 
of the instrument somewhat in order to maintain a suf- 
ficiently high sensitivity. 

The circuit of the level indicator is shown in Figure I (a). 
The sensitive element Rg is made of 0-15 mm diameter 
tantalum wire (resistance ~ 3-5 Q at 293°K and ~ 0-3 Q 
at 10°K) and is included in an arm of a bridge formed by 
resistances R; =2°6 Q, R,=14:5 Q, and the variable 
resistance R; = 2 Q. The supply from a 1-5 V source is fed 
through a button contact; a 10 mA milliameter is in the 
diagonal of the bridge. The bridge is balanced when the 
sensitive element is in vapour above the liquid; when it is 
immersed in liquid helium a great increase in current in the 
measuring diagonal takes place (up to 10 mA). 

If the indicator is removed from the liquid helium vessel 
with the button contact made, a current of over 25 mA 
flows in the reverse direction in the diagonal. A shunt 
resistance R= 0-55 Q is used to protect the measuring 
system from such possible damage; this reduces the useful 
signal insignificantly but reduces the reverse current to 
9mA. 

To achieve a still further decrease in the reverse current 
on lifting the probe out of the liquid helium vessel, an 
automatic reduction in bridge sensitivity was provided. 
A ballast resistance of nickel wire, Ry, in series with the 
voltage supply was placed on the probe above the sensitive 
element. When the probe is in the vessel of liquid the 
resistance of the nickel wire is small and the sensitivity of 
the indicator is at the maximum, but when the probe is 
lifted out the ballast resistance increases, producing a 
reduction in the voltage supply of the bridge. 


+ Received by PTE Editor 10 April 1960; Pribory i Tekhnika 
Eksperimenta No. 2, p. 190 (1961). 
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The tantalum sensitive resistance and the nickel ballast 
resistance are made in the form of 1:5 mm spirals and are 


Figure 1 


mounted at the end of the probe in channels of a Bakelite 
tube 4 mm in diameter and | mm thick (Figure 1(b)); the 
length of the probe is determined by the dimensions of the 
liquid helium vessel. 

The liquid helium level is determined by lowering the 
probe slowly into the vessel. The approach of the milli- 
ameter needle to the zero division shows that 40-50 mm 
remain before the surface of the liquid. The moment at 
which the sensitive element is immersed in liquid is fixed 
very clearly. 

The tantalum wire heats up rapidly to above 4:3°K on 
raising the probe out of the liquid, and after 1-2 sec the 
instrument is ready for a repeat measurement. 
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A 43 kgauss Superconducting Solenoid 


THE possibility of generating high magnetic fields by 
means of loss-free superconducting coils is of great techno- 
logical interest. For many years it was believed that 
because of the relatively low bulk critical fields of super- 


conductors superconducting solenoids would be restricted - 


to fields of a few thousand gauss. Recently, it has been 
realized' that by special treatment of materials, super- 
currents can be maintained at fields between one and two 
orders of magnitude greater than the bulk critical field. 
The first superconducting coil consciously to utilize this 
principle appears to have been wound by Yntema? using 
cold-worked niobium. This coil generated 7 kgauss. A 
similar coil was constructed by Autler,> who also estab- 
lished that superconducting windings could be used 
effectively in an iron-cored magnet. A further advance by 
Kunzler et al.* showed that 15 kgauss could be attained 
in a coil of cold-worked molybdenum-rhenium wire. : 

A great impetus to the development of superconducting 
coils occurred when Kunzler and his associates! demon- 
strated high supercurrent densities in short length Nb3Sn 
composite wires for externally applied fields up to 88 
kgauss. Despite the fact that this material is believed to 
carry supercurrents up to fields in excess of 200 kgauss, 
metallurgical difficulties in working with such brittle 
material in long lengths seem to have hindered the con- 
struction of actual solenoids. Since alloys would be ex- 
pected to be less brittle than intermetallic compounds, 
attention has turned to some of the high transition tem- 
perature alloys of niobium, in particular the niobium-— 
zirconium system.>’ Again, short length heavily cold- 
worked specimens of niobium-zirconium wire exhibited 
quite high supercurrent densities up to 30 kgauss.® 

Typical critical supercurrent densities for niobium— 
zirconium wire appear to be between 10* and 10° A/cm’. 
Thus, it can easily be calculated that to construct a small 
solenoid of 10 mil wire to generate a field of 50 kgauss, one 
needs several thousand feet of wire. This wire must be 
reasonably uniform in metallurgical characteristics since 
the current carrying capacity of the wire is limited by the 
lowest value encountered along its length. 

We have recently succeeded in producing niobium— 
zirconium wire of the desired uniformity in length up to 
several hundred feet. A 75 per cent niobium, 25 per cent 
zirconium base alloy was vacuum-arc melted and cast into 
4 in. diameter pins several inches long. These pins were 
canned in stainless steel and reduced by hot forging fol- 
lowed by cold swaging to } in. diameter rods. The rods 
were drawn to 0-010 in. diameter wire by the Armetco 
Company of Wooster, Ohio. 

Eighteen pieces of wire varying in length from 50 
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to 400 ft were insulated with several coats of epoxy 
enamel with a total build of 4 mil. The wire was wound 
on to a } in. diameter micarta phenolic resin former of 
length 1-4 in., resulting in a 74 layer coil of about 128 
turns per layer (Figure 1). Individual lengths of wire 


Figure 1. 43 kgauss superconducting coil 


were terminated at the end of a layer and brought 
out of the coil to a region of low magnetic field where a 
superconducting joint could readily be formed with the 
next piece of wire. The individual lengths were grouped 
into four sections consisting of 12, 24, 14, and 24 layers, 
respectively. Current was fed to each of the four sections 
from a separate individual power supply. 

Short test samples of wire prepared in the above manner 
were found to carry between 20 and 40 A at 20 kgauss and 
4-2° K. Individual samples frequently gave critical currents 
varying as much as 10 A between repeated tests. The cur- © 
rent carrying capacity of the four sections of the solenoid 
ranged from 15 to 20 A depending upon whether the 
sections were activated together or run up separately. We 
conclude from this that the wire is fairly uniform in long 
lengths. 

With the maximum permissible current in each section, 
a magnetic field slightly over 43,000 gauss was generated 
in the centre of the solenoid at 42° K. The field was mon- 
itored by a search coil and by a Hall probe at lower values. 
The field was found to be proportional to the current 
through the coil and in good agreement with the value 
calculated from the coil geometry. 

When the threshold current was attained in an individual 
section of the coil, the impedance rose very rapidly and 
reduced the current in that section to a fraction of an 
ampere. This frequently, although not always, caused the 
other sections to transit to the normal state, so that the 
entire peak field of 43,000 gauss was switched off in a 
fraction of a second. This behaviour subjected the coil toa 
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mechanical impulse which was clearly audible outside the 
metal cryostat in which the system was contained. Pre- 
sumably, this impulse is due to the sudden release of the 
forces acting upon the inner windings of the coil due to the 
high magnetic field. 


We are grateful to R. L. Ammon, W. Bagdade, J. H. 
Bechtold, R. T. Begley, R. D. Blaugher, G. Comenetz, 
R. E. Gold, G. T. Mallick, A. Patterson, S. W. Sand- 
berg, F. A. Sattler, and C. B. Satterthwaite for active 
assistance in parts of; this work. 
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EDITORIAL NOTE 


IN our first issue of September 1960 we mentioned that bibliographical material covering 
the low temperature field was being made available to our Journal through the kindness 
of Arthur D. Little Inc. of Cambridge, Massachusetts. In 1952 A. D. Little published a 
list of publications, each with a short abstract, comprising just under 1,000 entries, 
covering the years 1944-51. This material was broken down into a number of sub- 
sections and the publication was sent out to workers in the field of cryogenics and others. 
This first publication was augmented by nine supplements, giving titles of papers only, 
covering the years 1952-57. These supplements were in alphabetical order of authors. 
Although no further publication was made, A. D. Little kept their records up to date 
and kindly offered us the whole of the material. 

The publishers of CRYOGENICS accepted this offer and the whole bibliography has 
been re-classified and brought up to date by Dr. D. F. Brewer of the Clarendon Labor- 
atory. There are about 5,000 entries altogether, covering publications in the low 
temperature field from 1944 up to 1960 when the quarterly bibliographies published in 
CRYOGENICS began to appear. 

This bibliography has now been published as a supplement to CRYOGENICS. The 
entries, each giving the authors and title of paper, are arranged alphabetically under 
22 headings in five main groups. Copies are obtainable from: CRYOGENICS, Carlton 
House, Great Queen Street, London, W.C.2, price 25s. 
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